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3.1 INTRODUCTION 

3.1.1 Motivation for Transport Studies 

In the astrophysics community, ‘the solar flare problem’ 
is generally considered to be how to accumulate sufficient 
magnetic energy in one active region and to subsequently 
release it on a sufficiently short time scale. Satisfactory solu- 
tion of the solar flare problem will require at least two 
achievements by the solar physics community: first, con- 
vincing theoretical demonstration that one or more mecha- 
nisms of energy storage and release can occur; second, 
convincing observational demonstration that one (or more) 
of these theoretical processes actually does occur in the solar 
atmosphere. The contents of this Chapter essentially relate 
to the second problem, being largely concerned with how 
the energy released from magnetic form is transported 
through the solar atmosphere before escaping in the form 
of the radiant and mechanical energy signatures which we 
must interpret. 

A central point of general agreement concerning all 
mechanisms suggested for dissipation of magnetic energy in 
flares is that the actual sites of reconnection must involve 
scale lengths well below currently, or forseeably, achiev- 
able spatial resolution. Consequently, observational evidence 
in support of a flare theory is necessarily indirect, in the sense 
of not involving measurement of plasma parameters in the 
primary dissipation regions. Such indirect evidence may be 
of several kinds. Firstly, circumstantial evidence may be ob- 
tained by spatial resolution of the geometry, on a larger scale, 
of the magnetic environment in which the mechanism oper- 
ates. This may permit distinction between such options as 
emerging flux models and twisted arch models or between 
mechanisms driven by currents parallel to the magnetic field. 


as opposed to perpendicular. Spatial resolution also permits 
mapping of the paths of flare products. Secondly, temporal 
evidence can be obtained by use of high time resolution to 
set limits on instability growth rates, to imply the occurrence 
of repetitive or multiple dissipation, to indicate the produc- 
tion sequence of the various flare manifestations and, by 
causality arguments, to set an upper limit to the size of the 
primary dissipation site. Thirdly, thermodynamic evidence 
on the nature of primary dissipation is obtainable from the 
distribution, particularly in the impulsive phase, of the flare 
energy release over its various modes, i.e. fast particles, con- 
duction etc. For example the Petscheck mechanism releases 
a major fraction of the magnetic energy directly into bulk 
mass motion while the tearing mode initially results chiefly 
in plasma heating and particle acceleration. 

Transport of energy away from the primary sites, and 
its ultimate thermalization, depend not only on the primary 
mechanism itself but also on the larger-scale structure of the 
active-region atmosphere in which the transport occurs. Con- 
sequently the study of energy transport as a diagnostic of 
flare mechanisms involves extensive theoretical modeling of 
the transport processes, as well as observational input, to 
provide the framework for interpreting the observations. Use 
of such transport studies to infer properties of the initiating 
disturbance is essentially an inverse problem, and so carries 
the danger of indeterminacy through mathematical ill- 
posedness. For example, the thermal structure of a conduc- 
tively evolving atmosphere rapidly becomes only very weakly 
dependent on the heating function which initiated it, and so 
is a poor signature of this function. In such situations, the 
best strategy is to utilize jointly as many as possible inde- 
pendent signatures of the process to ndnirnize indeterminacy. 
It is just such a combination of independent signatures of 
the flare process (much more compelling when taken 
together) that the coordinated observational approach of 
SMM has rendered possible. 

In addition to these flare-oriented objectives, of course, 
the study of flare energy transport has contributions to make 
to the broader field of transport studies per se, such as in 
plasma and atomic physics. 

3.1.2 Historical Perspective 

Energy transport studies have become an increasingly 
prevalent means of investigating solar flares since around 
1970 with the accompanying steady improvement in observa- 
tional coverage and resolution in space, time, and spectrum, 
from the start of the Orbiting Solar Observatory (OSO) period 
onward. By the time of the Sky lab Apollo Telescope Mount 
(ATM), considerable progress had been made toward ob- 
taining the instrumentation needed for acquisition of high 
resolution data over as wide as possible a variety of 
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wavelengths from optical to 7 -rays. Thus the ATM package 
achieved spatial resolution of the order of arc seconds in the 
soft X-ray and ultraviolet ranges together with extensive UV 
line spectroscopy. Contemporarily though not simultane- 
ously, other satellites (notably ESRO TDl A, OSO-7 and the 
Intercosmos series) were improving the quality of hard X-ray 
measurements and extending the spectral range upward in 
energy resulting in the detection of solar 7 -ray lines. Though 
ATM itself was never designed as a flare mission it neverthe- 
less made major contributions to progress in flare observa- 
tions, in addition to its pioneering discoveries in relation to 
solar coronal structures and active regions (cf. Zirker 1977, 
Orrall 1981). In particular, ATM established the importance 
of loop structure in the magnetic configuration of many 
flares, and the compactness of bright XUV flare kernels. In 
addition, the use of ATM for flare studies (Sturrock 1980) 
delineated the limitations of such a package for answering 
some of the key questions concerning flares, and thereby 
provided important guidelines for the planning of subsequent 
missions, especially SMM. 

The most important obstacles to progress before the 
launch of SMM were the lack of data of sufficiently high 
time resolution, the lack of data sufficiently early in (or prior 
to) the flare, and the lack of data coordination over a wide 
enough spectral range (Brown and Smith, 1980). Examin- 
ing ATM in the light of these obstacles, we can see with hind- 
sight that ATM did not respond sufficiently quickly to enable 
systematic studies of flare onset or impulsive phases. Dur- 
ing ATM much information on flare spatial structure was 
recorded photographically in most characteristically ‘ther- 
mal’ wavebands — from optical to soft X-rays — with result- 
ing limitations on time resolution, on simultaneous spectral 
coverage, and on calibration. Typically ‘non-thermal’ emis- 
sions such as hard X-rays, 7 -rays, and microwaves were 
recorded with comparatively low sensitivity, and little or no 
coordination with the ATM experiments. A typical conse- 
quence of these problems was that testing of electron heated 
models of flare atmospheres involved construction of a spec- 
troscopic model atmosphere from data obtained from a va- 
riety of places in a variety of flares, and use of hard X-ray 
data, devoid of spatial information, from yet different flares. 

By contrast, SMM formed a coordinated package of in- 
struments dedicated to study impulsive phase phenomena at 
high time resolution using pre-planned targets and automated 
response to flare onset triggers. In addition, the spacecraft 
carried short wavelength instruments of unprecedented sen- 
sitivity and time resolution (the Gamma-Ray Spectrometer, 
GRS, Forrest et al. , 1980, and the Hard X-ray Burst Spec- 
trometer, HXRBs, Orwig et al.y 1980) and spatial resolu- 
tion (the Hard X-ray Imaging Spectrometer, HXIS, van Beek 
et al , 1980) and the facility for high resolution atomic X-ray 
line Spectroscopy (the Soft X-ray Polychromator, XRP, 
Acton et aly 1980), with digital data recording, as well as 
the Ultraviolet Spectrometer and Polarimeter (UVSP, 
Woodgate et al , 1980). Furthermore, the package was sup- 


ported by a wide range of other spaceborne instrumentation 
and an international network of ground based observations 
coordinated through the SMY (Svestka, Rust, and Dryer 
1982). Ground support included rapid arc-second resolution 
in both microwaves, by the VLA, and in spectrally resolved 
optical lines by the Sac Peak Vacuum Tower. Hard X- and 
7 -rays were still observed without spatial resolution but with 
such sensitivity as to permit close temporal correlation with 
features in the longer wavelength images, as well as 7 -ray 
nuclear abundance spectrometry with some simultaneous mil- 
limetric coverage at ultra high time resolution by Itapetinga. 
On the debit side, SMM was limited in soft X-ray spatial 
resolution compared to ATM, and in the restricted range of 
UV spectral coverage, of particular importance in modeling 
transition-region lines. 

A significant contribution to our knowledge of impulsive- 
phase transport has come about as a result of the post-SMM 
launch of the Japanese Hinotori spacecraft, whose flare in- 
struments (Solar Gamma-Ray Detector, SGR, Hard X-ray 
Monitors, HXM and FLM, Imaging X-ray Telescope, SXT, 
and Soft X-ray Crystal Spectrometer, SOX) are described 
by Kondo (1982). Early Hinotori results have been described 
in the Hinotori Symposium on Solar Flares (Tanaka et al , 
1982) and Recent Advances in the Understanding of Solar 
Flares (Kane et al , 1983); later Hinotori results play an im- 
portant role in this chapter. 

In addition to the observational requirements already men- 
tioned, there is obviously a need for improved theoretical 
modeling, particularly in the direction of making predictions 
which would be testable in terms of realizable data (Brown 
and Smith, 1980). The post- ATM period has indeed seen a 
major increase in the amount and sophistication of flare 
modeling work, particularly in respect of relaxation of earlier 
simplifying assumptions (such as hydrostatic equilibrium and 
optical thinness) in describing the atmospheric response, and 
of basic electrodynamic and plasma collective effects on the 
transport of charged beams. 

The perceptive reader will find that many of the ques- 
tions posed in the Skylab era have been answered in this chap- 
ter. Kane et al (1980), in the impulsive phase chapter of 
Solar Flares: A Monograph from Skylab Solar Workshop II 
(Sturrock, 1980), posed three key questions: 

1 . Is the distribution of energetic electrons thermal or 
nonthermal? 

2. Do the energetic particles (electrons), produced dur- 
ing the impulsive phase, provide the energy for the 
whole flare? 

3. Among the models of the impulsive phase suggested 
so far, which ones are most consistent with ob- 
servations? 

The answers to these questions are fundamental to our 
understanding of space plasmas; if the models that are best 
supported by the observations require particle acceleration 
efficiencies >0.1% (Hudson, 1979), a substantial challenge 
is presented to the solar flare theorist. 
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Our work answers primarily the first and third questions, 
while the second motivated the Solar Flare Energetics group. 
Certainly there are still major gaps in our theoretical under- 
standing of how energy is propagated by electron beams, as 
well as our observational understanding of spatial, temporal 
and spectral scales. However, one impulsive phase transport 
model now stands out above all others: the nonthermaU 
electron thick-target model, in which the dominant role in 
the transport of energy on impulsive-phase timescales (usu- 
ally ^ tens of seconds) is played by beams of electrons, 
mostly in the deka-keV range, whose velocity distribution 
function cannot be described by a single-temperature Max- 
wellian. These electrons are guided along a loop-like mag- 
netic field structure, from an acceleration site in the corona; 
they heat and cause both thermal and nonthermal emission 
as they are fully thermalized in the loop plasma, a thick tar- 
get. The reader should not get the impression that this model 
describes all flare energy transport, particularly on longer 
timescales, or that the nonthermal-electron thick-target model 
passes all the observational or theoretical tests we impose 
upon it below. Surely there is ample evidence for the 
impulsive-phase existence of nonthermal particles other than 
electrons, for example, and for thermal domination of later 
flare phases. However, our preoccupation with the 
nonthermal-electron thick-target model in this chapter is testi- 
mony to our finding that, of the models available, it does 
the best job of explaining the wide variety of impulsive-phase 
observations we have studied. 

3.1.3 Overview of the Chapter 

Some of the most striking recent results have come from 
observations in the highest energy ranges. Interpretation of 
SMM 7 -ray line data has established the presence of pro- 
tons (E ^ 10 MeV) in regions of density ^ 10^^ cm”^, 
while the ISEE-3/PVO occultation data demonstrated that 
electrons of E ^ 150 keV are stopped deep in the chro- 
mosphere. Both these results suggest a thick target beam in- 
terpretation, in which the radiation is generated in the course 
of fully stopping the particles, in the lower solar atmosphere. 
On the other hand, ISEE-3/PVO data surprisingly show very 
little directivity at 350 keV, contrary to purely collisional 
transport models of the electron beam. Most notable of all 
high energy data is the striking demonstration, from limb 
brightening studies, that directivity is present in the con- 
tinuum around 10 MeV. 

The first-ever images in the deka-keV range (10 - 100 
keV) stimulate much of the work of this Chapter. While there 
remains considerable debate over their implications for theo- 
retical models (in particular beams), the presence in some 
flares of impulsive hard X-ray footpoints, coincident in time 
and space with chromospheric emissions, has been clearly 
established by SMM and Hinotori. Hinotori results have also 
been used to suggest a provisional classification of hard X-ray 
flares, of which footpoint events are only one class. Simul- 


taneous microwave data have permitted comparative mor- 
phology studies of hard X-rays and microwaves with 
resulting constraints on the electron and magnetic field dis- 
tributions. Computations of hard X-ray polarization incor- 
porating the effects of magnetic field curvature show that 
even the most recent low polarization results from Shuttle 
experiments are not incompatible with a beam model. 

While the greatly improved spatial and temporal resolu- 
tion and coordination of SMM data in the hard X-ray and 
EUV further support the view that these radiations have 
closely related origins, the relationship between the fluxes 
in these bands remains a theoretical enigma in terms of 
energy balance. 

The inclusion of high resolution X-ray spectroscopy al- 
lowed detailed diagnostics of the hot plasma during the on- 
set of the thermal phase, using recent developments in the 
atomic physics of dielectronic satellite spectra. Line profiles 
and shifts enabled the determination of turbulent energy con- 
tent and upward velocities during the in^ulsive phase, result- 
ing from the chromospheric evaporation process. 
Observations of X-ray and Ka line excitation, together with 
data on the soft X-ray plasma as a whole and on the hard 
X-ray source electrons, have demonstrated that electron beam 
excitation is not essential to explain the Ka observations, 
though not precluded by it. 

Improved calculations of the production of the non- 
thermal red shifted La line by capture processes on a 
descending proton beam have placed severe constraints on 
the flux of such beams and hence on flare models where pro- 
tons play a central role energetically and in which there has 
been renewed interest recently. 

At the time of ATM, white light flares were stiU regarded 
as a rare phenomenon suggestive of an exotic explanation. 
Studies with improved observational methods have shown 
that they are in fact of common occurrence. Coordinated tem- 
poral, spatial, and spectral information shows that the im- 
pulsive phase white light emission is spectrally compatible 
with a chromospheric origin, and temporally and ener- 
getically compatible with a thick-target electron beam as the 
power supply. 

Interpretation of Ha profiles represents a particularly 
good example of the progress SMM achieved toward the 
ideals of rapid response, data coordination, and detailed 
modeling. Full radiative transfer models of the flaring chro- 
mosphere have permitted comparison on a pixel by pixel, 
instant by instant basis of Ha profiles observed at Sacramento 
Peak widi those predicted according to chromospheric evapo- 
ration models driven respectively by thick target electron 
beams (observed by HXIS) and by thermal conduction. 
Results provide direct evidence in support of electron beam 
heating in some bright Ha kernels. 

Finally, prior to the SMM era decimetric radio waves 
were considered to be a plasma diagnostic emission of little 
importance energetically. At this time serious consideration 
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is being given to their role in heating both the corona (by 
reabsorption) and the chromosphere (by electron precipita- 
tion) in flares where conditions favor the onset of coherent 
amplification of these radio waves. The full ramifications 
of this process for flare energy transport are only now be- 
ginning to be appreciated. 

3.2 IMPULSIVE PHASE OBSERVATIONS 
AND THEIR INTERPRETATION 

3.2.1 Gamma-Ray Emission Above 10 MeV. 

Photon emission above 10 MeV during solar flares was 
observed by the Gamma Ray Spectrometer (GRS) on board 
SMM. As of February 1984 the highest energy photons de- 
tected, energy - 80 MeV, came from an intense flare on 
June 3, 1982. In this contribution we discuss the timing be- 
tween hard X-rays and gamma rays and present evidence for 
directivity of the highly energetic particles that give rise to 
the emission above 10 MeV. 

3.2.1. 1 Relative Hard X-ray and Gamma-Ray 
Timing 

As a typical example we show in Figure 3. 1 the time his- 
tory of the flare of June 15, 1982 in different energy bands. 
The event has a simple time structure consisting mainly of 
one impulsive burst. The peak of the emission is simultane- 
ous within ± 2 sec over more than three decades of energy 
from 30 keV ~ 50 MeV. The hard X-ray flux shown in the 
upper two panels is assumed to originate from electron 
bremsstrahlung. The emission from 4.1 - 6.4 MeV (the 
nuclear energy band) is from nuclear lines and from brems- 
strahlung of relativistic electrons (Forrest 1983). At ener- 
gies greater than 10 MeV the gamma rays are expected to 
be produced by bremsstrahlung of very highly energetic elec- 
trons and by the decay of pions (tt® and x"^), because the 
contribution of nuclear lines above 8 MeV is negligible 
(Crannell, Crannell, and Ramaty, 1979). The relative im- 
portance of the two processes for solar flares is discussed 
by Ramaty et al. (1983) and by Rieger et al (1983). 

Until February 1984, 14 flares were observed with emis- 
sion above 10 MeV. Their time history has the following 
general characteristics: They are of short duration ( — 1 min) 
and very impulsive. Rise and fall times are on the order of 
seconds. They exhibit single or multiple peaks. The peak 
emissions of the hard X-ray and gamma rays are simultane- 
ous within about ± 6 sec. There is no systematic delay of 
the gamma rays with respect to the X-rays. 

Because of the simultaneity of the gamma- and X-rays, 
the energy loss of the highly energetic particles (electrons 
and/or ions) has to take place at ambient densities > 10*^ 
cm~^, if we assume a thick target situation (Bai and Ramaty 
1976). 
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Figure 3.1 Time history of the flare of June 15, 1982 
in various energy bands. 

3.2. 1.2 Directivity of Highly Energetic Particles 

The flares with photon emission above 10 MeV are at 
heliocentric angles of > 60° (Rieger et al, 1983). This is 
shown in Figure 3.2, where the location of the flares, known 
from Ha observations (NOAA: Solar Geophysical Data) is 
plotted. The arithmetic mean heliocentric angle of all 14 
flares is 79 °. If the radiation is isotropic the probability for 
a chance coincidence of such a distribution is ^ 2 X lO""^. 
Therefore we conclude, that this “limb brightening” is the 
result of directivity of the radiation. This directivity must 
exist also in the primary highly energetic particles, because 
photons with energies above 10 MeV, if created by electrons 
( > 10 MeV) via bremsstrahlung or by protons ( > 100 MeV) 
via pion decay, are emitted preferentially in the direction of 
the motion of these particles. At this time the most likely 
interpretation appears to be that the emitting particles are 
indeed travelling roughly normal to the sun-center direction 
(J. Cooper and V. Petrosian, private communication), as if 
they were near their mirroring point, for example. 

To study this phenomenon in more detail it would be 
necessary to make stereoscopic observations with two de- 
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Figure 3.2 Solar disc position of flares with photon 
emission above 10 MeV. 

tectors of the GRE type widely spaced in solar aspect angle, 
as was done by Kane et al (1980, 1982) at lower energies. 

3.2.2 Hard X-ray and Microwave 
Morphology 

3.2.2. 1 Observations 

The strong correlation in the temporal evolution of hard 
X-rays and microwave radiation established over the years 
has led to the belief that the same or very closely related 
populations of electrons are responsible for both of these two 
different flare radiations. High spatial resolution observa- 
tions at both wavelength ranges provides bjrther informa- 
tion on the details of this correlation. 

Since the first fanbeam observation of flares by Enome 
et al. (1969) with resolution 24" at 9.4 GHz, the new in- 
struments such as the NRAO 3 element interferometer with 
resolution 10" at 2.7 GHz (Alissandrakis and Kundu 1975), 
the WSRT array with one dimensional resolution of 6" at 
5 GHz (Alissandrakis and Kundu 1978, Kattenberg and 
Allaart 1981) and finally the VLA with two dimensional reso- 
lution of 1 " at 15 GHz (Marsh et al. , 1980, Lang et al. , 1981 
and Kundu et al., 1981) have steadily improved the resolu- 
tion of microwave observations. VLA observations (Kundu 
et al., 1982, Dulk et al., 1983) and fanbeam observations 
at 35 GHz (Nagoya) and 17 GHz (Nobeyama) have con- 
tinued, while the HXIS on SMM and the SXT on Hinotori 
have begun for the first time to provide spatially resolved 
images with resolution of about 10 " at X-ray energies of less 


than 40 keV. Unfortunately, simultaneous microwave and 
hard X-ray high resolution images exist for only a few flares. 
Consequently, to some extent, we must still rely on the data 
of more numerous events with only either hard X-ray or 
microwave images. 

We first summarize the results of observations where 
there exists only high resolution microwave data, then those 
with only X-ray data and, finally, the few cases with both 
X-ray and microwave observations. 

a. Microwave Morphology. Studies of spatial structure 
of microwave radiation from flares have been summarized 
by Marsh and Hurford (1982) and Kundu (1983). It should 
be noted at the outset that, unlike the X-ray studies, the 
microwave observations (in particular the high resolution 
one-dimensional results from WSRT and the two-dimensional 
results from the VLA) have not been as systematic, continu- 
ous and extensive as the X-ray observations. Consequently, 
it is difficult to classify the microwave structure in well de- 
fined categories. In spite of this, these observations have 
shown some common features on which we will concentrate 
here, remembering that there may be more complex struc- 
tures yet to be studied, analyzed, and classified. 

The most general statement that can be made is that the 
brightest point of the microwave radiation during the impul- 
sive phase occurs near a magnetic neutral line and not on 
an Ha kernel, and that in cases with simple field geometries 
there is a single dominant source. Whenever a secondary 
source is detected, that source lies near another neutral line 
(e.g., 26 May 1980 flare, Kundu 1983). It can be concluded 
that with few exceptions (see Kundu et ai, 1982), the 
predominant microwave emission does not come from the 
footpoints of the flaring loop. In some flares the microwave 
source is approximately midway between the footpoints and 
is nearly equally but oppositely polarized (circular) on both 
sides of the maximum brightness spot (see Figure 3.3). The 
simplest interpretation of these observations is that the emis- 
sion comes from the region around the ton of a loop. Spec- 
tral variation of the structure and polarization then tell us 
something about the geometry of a loop and the pitch angle 
distribution of the radiating electrons (Petrosian 1982). 
However, there are flares where only one sense of polariza- 
tion is observed (cf. Kundu 1983), and there are cases where 
the source is not located very high up in the corona (Kai et 
al. , 1982). These features have b^n attributed to asymmetric 
magnetic loops (Kundu and Vlahos 1979) but, as we shall 
discuss below, such structures could also arise when a non- 
circular or sheared loop is viewed from an angle away from 
the line perpendicular to the field at the top of the loop. 

b. X-ray Morphology. Prior to the era of SMM and 
Hinotori, the only information on hard X-ray spatial struc- 
ture was obtained from stereoscopic observation by the PVO 
and ISEE-3 satellites (Kane et al., 1979, 1982). These ob- 
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Figure 3.3 VLA map at 15 GHz in right circular (black) and left circular (white) polarization superimposed on 
Ha frames. The dashed line depicts the location of the magnetic neutral line (from Hoyng et a!., 1983). 

servations do not provide images, but give easily-interpreted There are also rarer occurrences of double sources (not 

information on the height of the hard X-ray emitting regions. necessarily of equal intensity) which merge into a single 

Furthermore, they extend to much higher X-ray energies than source located between the original two sources as the flare 

the imaging instruments and tell something about the varia- progresses. Similar structures and evolution have now been 

tion with height of X-ray spectra, which turns out to be very shown to be present in the HXIS data (Machado 1983. and 

useful in modeling (cf.. Brown et. al., 1981, Leach and Machado el al, 1983). 

Petrosian 1983). Based on light curves and spectral images measured by 

The first hard X-ray images of flares were obtained by instruments on board Hinotori, various classifications of the 

the HXIS on SMM. The analysis of the few early flares em- flares have been given by Ohki et al. (1983). Tanaka (1983) 

phasized X-ray structure consisting of two sources which and Tsuneta (1983). For the following rea.sons. the classifi- 

were identified as the footpoints of flaring loops (Hoyng et. cations are preliminary at this time. They are based on only 

al., 1981). But as indicated more recently (Duijveman and the strongest flares, which may not be representative of all 

Hoyng 1983), there are many flares where the bulk of the flares. The majority are weaker and shorter-lived than the 

hard X-rays (> 16 keV) come from a single .source. ones used for this classification. In addition, relatively few 

Analysis of the data from SXT on Hinotori (Ohki et al. , flares have been classified; the number of flares in each cla.ss 

1982, Tsuneta et al . , 1982, Takakura et al. , 1982, Ohki et is small. It should also be noted that they are based on only 

al., 1983, Tsuneta e/ «/., 1983, Takakura e/ a/. , 1983 and qualitative differences between spatial structures, spectra and 

Kosugi etal. , 1983) indicates the predominance of flare im- the impulsiveness of the light curve. For a more thorough 

ages consisting of a dominant .single source which sometimes analysis the parameters describing these characteristics 

expands and become elongated, and sometimes shrinks. should be quantified and flares binned accordingly. The flares 
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studied may be the extreme cases, so that further analysis 
may show a continuum of classes with no clear dividing line 
between the present categories. 

The Hinotori data have been classified into three types, 
A, B, and C, as described in Table 3.1. Type A events are 
defined as those which have a smooth light curve below 40 
keV (no discernable spiky features), a soft spectrum (spec- 
tral index > 6 for a power law fit, or an exponential spec- 
trum within observational uncertainty) and consist of a single 
compact source. Two examples of such flares have been ex- 
tensively studied (Tsuneta et aL, 1984), the April 2nd and 
July 17th flares of 1981. The April 2nd flare (Figure 3.4) 
shows a single steady point source (perhaps barely resolved 
at 15 ")• Strong Fe XXVI lines are also observed throughout 
these events. 

Type B flares are those with two distinct hard X-ray bright 
spots during initial (impulsive) part of the flare, which evolve 
into a single source located somewhere between the initial 
bright regions. During this evolution it is observed that spiki- 
ness of the light curve disappears and the spectrum softens. 
Hinotori events of July 20 (Figure 3.5) and October 15 of 
1981 are two such events. 

Finally, there are Type C flares which, like Type A flares, 
consist of a single source, but the source is more diffuse and 
is clearly displaced from the Ha kernel and, like the micro- 
wave sources described above, could be emanating from the 
top of a set of loops. The light curve is also smooth, even 
at energies greater than 49 keV, and the spectrum is hard. 


The May 13 (Figure 3.6) and April 27 flares of 1981 are 
two such events. 

As mentioned above, the HXIS data also show flares with 
spatial structures similar to those described by Hinotori. 
Figure 3.7 shows some such examples. 

It should be noted that some of the Type A flares may 
actually be of Type B, but unresolved, either because they 
are intrinsically compact or because they are elongated 
sources viewed along the major axis. 

c. Simultaneous Hard X-ray and Microwave Morphol- 
ogy. There are less than a dozen flares which are resolved 
(in one-dimension or two) at both hard X-ray and microwave 
energies. Here we describe the brightest of these and con- 
centrate on the relation between the hard X-ray and micro- 
wave structure of the most important features. We begin by 
summarizing Hinotori observations of four events represent- 
ing the three types described above, and then four SMM 
X-ray events that have also been observed in microwaves 
by the WSRT or the VLA, and fit them into categories A, 
B or C. 

The August 1 1 , 198 1 , Hinotori flare is classified as Type 
A, but as far as the X-ray image is concerned, it falls on 
the border line between Types A and C. Without a quantita- 
tive measure of the distinguishing parameters of the various 
types, one cannot resolve this ambiguity. In any event, the 
X-ray image appears to be a steady single elongated source 
(with a hint of a double source at the rising phase). The one- 


Table 3.1 Hinotori Flare Classification 

OBSERVATIONAL CHARACTERISTICS 
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Figure 3.4 Example of a Type A flare, Aug, 11,1 981 . 
The lower right-hand panel shows a sketch of the Ha 
flare at 01 :55:45 UT and the direction of the peak of 
the one-dimensional brightness at 17 GHz at 
01 :45: 1 0 UT. The other five panels show five differ- 
ent X-ray (25-50 keV) images (Takakura etaL, 1983). 

dimensional 17 GHz tracing shows a near coincidence of the 
peak microwave and X-ray positions (Figure 3,4). Both max- 
ima seem to be located over a neutral magnetic line. (For 
details see Takakura et ai, 1983.) 

The April 1 and May 13, 1981, Hinotori flares are both 
classified as Type C. The May 13 flare has a very smooth 
time profile even above 1(X3 keV, while the X-ray light curve 
for the April 1 flare is midway between the August 1 1 and 
May 13 flares. As shown in Figures 3,6 and 3.8, the 35 GHz 
microwave position agrees with the position of a single 
dominant X-ray source. Both emissions appear to come from 


the top of a loop or arcade of loops. However, there is some 
asymmetry in the one-dimensional microwave image. This 
asymmetry becomes stronger (indicating the existence of two 
sources) in the later phases of the flare, and in the case of 
the April 1 event (Figure 3.8) it may be interpreted as a 
double source. Unfortunately, there is no X-ray image for 
these periods. 

The October 15, 1981, Hinotori flare may be classified 
as Type B because initially it consists of two sources (one 
much brighter than the other). However it is not clear if these 
sources correspond to any Ha kernels. The 35 GHz image 
can be decomposed into a double source with the stronger 
microwave source corresponding with the weaker X-ray 
source (Figure 3.9). There is a secondary peak in the light 
curve during which the X-ray structure changes rapidly. The 
brighter source moves about 30 arcsec southward and the 
weak source becomes stronger. Finally, during the decay- 
ing phase these sources disappear and a third source appears 
in a new position. There are no microwave images for the 
second peak or the decay phase of this flare. Probably this 
is a Type B flare like the July 20, 1981 , limb flare. However, 
one cannot rule out the less likely possibility that the four 
different bright regions correspond to the top of four differ- 
ent loops. 

The November 5, 1980 SMM event (Figure 3.10) 
described by Hoyng et al. (1983) has excellent overlapping 
HXIS and VLA coverage. As is typical, the 15 GHz radia- 
tion comes from above a neutral line, as shown in Figure 
3.10, The hard X-ray and 9.4 GHz light curves show a Type 
B time profile and consist of three peaks, with the third peak 
being a soft, gradual one. The 16 to 30 keV HXIS counts 
in various pixels conform to the Type B character of this 
event in that, initially, there is less flux from the region over 
the neutral line than from the Ha bright spots, and during 
the second and third peaks most of the X-rays come from 
the approximate center of gravity of the bright region and 
coincide with the bright microwave source. 

The July 13, 1980, SMM event described by Kattenberg 
et ai (1984) shows two sources in coincidence with the Ha 
bright patches. Unlike the Type B flares there is no gradual 
peak (or phase) in the light curve and no emergence of a 
single source in the middle. Unfortunately, the baseline of 
the WSRT one-dimensional microwave tracing is perpendic- 
ular to the line connecting the two X-ray sources so that we 
do not know if the microwave source is also double. 
However, the usual centrally-located microwave image is 
consistent with the data. 

On June 24, 1980, two bursts at 15:20 and 19:57 were 
extensively mapped at 6 cm by the VLA with some cover- 
age by the HXIS on SMM (Kundu et al . , 1983). The 15:20 
event shows a hard X-ray image displaced by 20 to 30 arc- 
sec from the microwave source, the latter being typically 
bipolar (separated right and left circularly polarized struc- 
tures, as in Figure 3.3). The softer X-ray and microwave 
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UT 13:15:45 - 13:15:51 


UT 13:24:48 - 13:24:55 



Figure 3.5 Example of a Type B flare, July 20, 1981: (a) impulsive phase, (b) gradual phase. Each image is 
3' In size and each pixel is 6". The curved line bisecting the images shows the location of the solar limb. Note 
the evolution from a lower double to a higher single source (from Tsuneta et aL, 1983b). 


emissions are co-spatial, presumably coming from the loop 
top, while the hard X-rays may come from a footpoint. This 
may therefore be a Type B event. 

d. Summary. The sample of flares described above is ob- 
viously incomplete in many ways. In general, any sample 
is limited by selection effects associated with observing in- 
struments. However, certain inferences about flare phenom- 
ena are affected less than others by such selection criteria. 
We hope that by concentrating on the gross features of the 
flare phenomenon the influence of the selection effects can 
be minimized. 

The prominent features of the X-ray and microwave mor- 
phology on the basis of the small number of events presently 
analyzed, can be summarized as follows: 

1 . For almost all flares where there exist simultaneous 
microwave and Ho images the microwave emission 
comes from a region above the neutral line (presum- 
ably around the top of a loop) between Ha bright 
patches (presumed to be the footpoints of loops). 
There are exceptions, when a minor microwave fea- 
ture may be near an Ha bright patch (Dulk et aL , 
1983). Also, there are flares (e.g., Oct. 10, 1981) 
where because of coincidence of the microwave and 
hard X-ray images, it is assumed that all the emis- 
sion is from footpoints. However, there is no direct 
evidence for microwave emission from the footpoint 
of loops. 


2. In the majority of cases the X-ray emission comes 
from a single dominant source (Types A and C), 
which is found to be located above a neutral line and 
not on an Ha bright patch, except that for some com- 
pact, Type A bursts higher resolution is needed for 
confirmation of this picture. The structure of Type 
B bursts with two bri^t X-ray sources (normally coin- 
cident with Ha patches) evolves into a single source 
located between the Ha patches. 

3. The limited hard X-ray and microwave data obtained 
so far are consistent with, but are not a direct evidence 
for, the conclusion that such large scale changes in 
the images occur on a hydrodynamic time scale, with 
velocity v « 50 to 200 km/sec. 

4. The hard X-ray and microwave images, even though 
different in detail, have roughly the same location 
whenever both images are dominated by a single 
source. For bursts whose X-ray images are initially 
double, the single source that develops later most 
probably coincides with the microwave source, 
presumably located near the top of a loop or arcade 
of loops. 

As we shall show in the next section, the above general 
features are consistent with the nonthermal model, whereby 
semi-relativistic electrons are injected in a closed loop some- 
where in the corona. One aim of such observational and theo- 
retical work is to determine what constraints the observations 
provide on model parameters such as field geometry, loop 
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Figure 3.6 Example of a Type C flare, May 13, 1981. 
Comparison of the hard X-rays (17 to 40 keV), the 
Ha photograph and the one-dimensional radio images 
for the Indicated times (from Tsuneta et at., 1 983b). 


size and density, and the spectrum and pitch angle distribu- 
tion of the accelerated electrons. 

3.2.2.2 Interpretation 

We now interpret the above data in the framework of the 
nonthermal models whereby electrons are accelerated to ener- 
gies well beyond the thermal energy of the coronal plasma. 
These models are variously classified as thick target, thin 
target, trap, etc. 

The traditional trap models require both a rapid conver- 
gence of the field lines and a low density. For electrons of 
energy E to mirror back and forth in the corona, one needs 
6B/B > > 6L/X, where 6L is the length scale of the mag- 
netic field variation 6B, and X is the mean free path of elec- 
trons of energy E. For the X-ray -producing electrons 


(energies E — 20 keV) to survive for time t one needs a value 
of density n < W (E/20 keV)3/2(ioo/t sec) cm-\ 

The thin target model will be realized only in an open 
magnetic field configuration. Isotropic or outward stream- 
ing injections of electrons can be dismissed because, con- 
trary to observations (Lin and Hudson 1971, 1976), they 
require as many electrons to reach the earth as are needed 
to produce the hard X-rays. If the electrons are highly beamed 
toward the chromosphere, the result will be similar to the 
thick-target closed-field configurations discussed below. 

Consequently, we will consider only the general thick- 
target model where electrons are injected in a closed mag- 
netic loop (Figure 3,11), We shall need to specify the site 
of injection of the electrons. If the site of injection is very 
deep, well below the transition zone (point 1 of Figure 
3. 1 la), the electrons lose energy quickly, causing direct heat- 
ing and evaporation of chromospheric plasma, and produce 
hardly any hard X-rays above the transition region. We shall 
not discuss this possibility here because it will fail to describe 
Type C bursts discussed above. On the other hand, if the 
injection is at lower densities (point 2 of Figure 3.11a), so 
that all of the outgoing particles eventually reach the top of 
the loop, then the situation is qualitatively similar to the in- 
jection at the top of the loop (point 3), which is what we 
shall assume. 

a. Transport of Electrons, The transport of the nonther- 
mal electrons is affected by many parameters. In particular, 
if the pitch angle distribution of the electrons is highly 
anisotropic (beamed electrons), they are subject to some in- 
stabilities which change the distribution on a time scale of 
the order of plasma oscillations. We shall not consider such 
cases here. Furthermore, if the beams constitute a high cur- 
rent, then a reverse current is set up. The electric field which 
drives the reverse current may then decelerate the original 
beam. These aspects of the transport are discussed in Section 
3.3.1. Here we concentrate on the collisional effects and the 
role of large-scale static magnetic fields. 

The solar flare plasma conditions allow for various sim- 
plifying assumptions, so that the variation with depth of the 
distribution of electrons F(E,/i,r), where F dE d^i dr is the 
electron flux in the energy range dE (in units of mc^), pitch 
angle cosine range dfi and dimensionless column depth range 
dr = dN/No (Nq = 5 X lO^^cm”^ is the column depth re- 
quired for stopping an electron with E = 1) is determined 
by a single parameter dfnB/dr, which describes the varia- 
tion of the magnetic field B along the loop, with respect to 
the distribution of plasma. Parameters of the injected elec- 
tron spectrum F(E,/a,t) = are the spectral 

index 6 and the pitch angle dispersion a©- Leach and Petrosian 
(1981) applied the Fokker-Planck method to electron trans- 
port in flare loops. Here we point out an aspect of their work 
that is relevant to the relationship between hard X-ray and 
microwave emission. 
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9:18:38 


JULY 14. 1980 
8:26:02 


MAY 9. 1980 



Figure 3.7 Examples from HXIS of double sources, single sources and evolving structures similar to those seen 
by Hinotori (from Duijveman and Hoyng 1983 and Machado 1983). 


Let us initially consider loops with uniform magnetic 
field, i.e. dfnB/dr = 0. For nonrelativistic electrons, pitch 
angle diffusion and energy loss take place on comparable 
time scales. This means that models with widely different 
values of aj appear different only at the top of the loop. Such 
models are virtually indistinguishable models lower down, 
where the electron distributions become nearly isotropic. 
Consequently, the radiation signature of spatially unresolved 
loops will tend to be closer to what one expects from an 
isotropic rather than a highly beamed pitch angle distribu- 
tion. Note, however, that this is not true at relativistic ener- 
gies (i.e. electrons producing the high-frequency microwave 
radiation and the continuum y-rays). Relativistic electrons 
lose energy more quickly than they diffuse in pitch angle. 

The effect of magnetic field convergence, i.e., of a non- 
zero dfnB/dr, is to mirror some of the electrons which do 
not reach to the chromosphere and confine them to the top 
of the loop. However, such electrons also eventually scatter 
into smaller pitch angles, penetrate deeper, and thermalize. 
The overall effect of the non-zero dfnB/dr is to isotropize 
the pitch angle distributions more quickly and enhance the 
effects described in the previous paragraph. 

b. Non-thermal Emission. The most prominent aspects 
of the hard X-ray and microwave observations can readily 


be understood in the framework of a simple semi-circular 
loop model, as described above. Although this model is too 
simplified to use in interpreting the observations in detail, 
it shows that the geometry of the loop and the physical 
parameters of the plasma play a more significant role than 
emphasized in previous treatments. Following Leach and 
Petrosian (1983) and Petrosian (1982), we assume a semi- 
circular loop (radius R) in the corona with variable magnetic 


mosphere (see Figure 3.11a). We also assume that the den- 
sity of the preflare plasma in the loop is constant in the 
corona and rapidly increasing in the transition region. The 
important parameters for our consideration are the column 
depth of the transition region, = noxR/2, and dfnB/dr. 

We begin with the spatial structure of X-rays in the range 
16 to 50 keV. The X-ray emission is obtained from the elec- 
tron flux F given by the transport analysis, and the brems- 
strahlung radiation cross-section. At these energies the 
electrons are non-relativistic, and their radiation is nearly 
isotropic. The scaling between electron energy E and colunm 
depth T directly translates into a similar scaling between r 
and the photon energy k (in units of m^c^). Leach (1984) 
showed that the variation of the X-ray intensity I (normal- 
ized photon counts per energy interval dk and depth interval 
dr) with k and r can be described to within 20% by 
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(IN SOLAR OPTICAL RADIUS) 

Figure 3.8 One-dimensional microwave tracing and 
the hard X-ray image from SXT of the April 1 , 1 981 , 
flare (from Kawabata et a/., 1982). 


I(k,T)dr « (l-hT/k2)5/2 dr/k2. (3.1) 

This simple picture becomes complicated when we trans- 
late it into observables. Let us consider the circular loop of 
Figure 3.11a viewed from directly above. For comparison 
with observation we need the variation of the intensity 
projected on the solar disk I(k,r) with the projected distance 
r from the center of the loop: 


15 OCT 1981 

\ 

SXT X-RAY IMAGE 



W0.30RO W0.40RO 


W5' W6' W7' 

Figure 3.9 A comparison of the SXT hard X-ray image 
with the fanbeam radio brightness after restoration for 
the October 15, 1981, event (from Kawabata et aL, 
1 982). Whether the two images belong to footpoints 
of a single loop or emanate from top of the distinct 
loop is not known. 



Tb l = 5.4x1O0 K 
VLA 15 GHz Tb*r = 9.4x10®K 


I(k,r) = I(k,r)dr/dr = I(k,r) (n^/N^) (l-rVR 2 )-i /2 (3 2 ) 

The last geometric term plays a significant role, as shown 
by ds/dr in Figure 3.12. Note, however, that even for non- 
circular geometry ds/dr (and dr/dr) can increase rapidly as 
the footpoints are approached. At a given photon energy k, 
if Njj. is small, or more precisely if < < k^, 

then I(k,r) will be nearly constant for 0 < r < R. However, 
because of projection effects I(k,r) will rise with r, such that 
the footpoints at r = R will appear brighter than the top of 
the loop at r = 0 (cf., Figure 3.12b). For > > k^N^ 



HXI5 16-30 KeV 


Figure 3.10 The microwave and hard X-ray images 
of November 5, 1980, flare. The lower panels are 
8" X 16" in size and a dot indicates significant flux 
in excess of a single temperature fit (from Duijveman 
and Hoyng 1983). (cf. Figure 3.3) 
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Figure 3.11 Description of the geometry of loops: 
(a) semi-circular loop; (b) asymmetric sheared loop. 
Circled points (1 ), (2) and (3) refer to the possible lo- 
cations of injection of electrons. The wavy arrows 
show the line of sight and the value of 0 used in equa- 
tion (3). Note that ds and di = n ds/No are measured 
along the loop from the point of injection, assumed 
to be at (3). R is radius of the loop and r is a projected 
distance from the center of the loop. 



I 



Figure 3.12 (a) Variation with projected distance (half 
of the loop) from the loop centers of various quanti- 
ties. (b) Variation of expected hard X-ray brightness 
for the indicated values of N^^/Nok^. 


the brightness will decrease steadily from the top to the foot- 
points. The changeover of the location of the maximum 
brightness from r = 0 to r = R will occur quickly around 
= 1Q20 cm“2(k/22 keV)^. 

At photon energies less than 10 keV(k < .02) a low pre- 
flarc value of N^j. — lO-'cm”* is sumcicnuy high to ensure 
that the maximum brightness occurs at the top of the loop. 
But such a low value of is not sufficient to confine 
higher-energy (say > 22 keV) electrons to the top of the 
loop. Such electrons wUl penetrate below the transition region 
so that their radiation at 22 keV will be concentrated on the 
footpoints. Only if > lO^^cm'^ does the top becomes 
brighter than the footpoints at these energies also. 

This general picture is nearly independent of the view- 
ing angle and the field geometry because in general dr/dr 
becomes very large near the footpoints. However, the de- 
tails of the shapes of the curves drawn on Figure 3.12 do 
depend on the field geometry and the viewing angle. 

This behavior therefore suggests that for the large diffiise 
Type C flares described in the first part of this section, 
N^j. > 10^0 cm~2 (perhaps n - 10*® cm"^ L -- 10*® cm), 
so that one sees the 10 - 40 keV photons from the tops of 


the loops. On the other hand, for the Type B events, ini- 
tially Njj < 10^® cm”^ and the emission is concentrated in 
the footpoints. The asymmetry of the radiation from the two 
footpoints as commonly observed must then be attributed to 
an asymmetric field geometry (e.g., Figure 3.11b). As ob- 
served by Hinoiori (Figure 3.5) and in some events by HXIS 
(Figure 3.7), the double images merge into a single coronal 
source as the flare progresses. This means that the column 
depth has increased from its initial low value to a value larger 
than 10^® cm~2. As mentioned in connection with observa- 
tions, such changes seem to occur on a hydrodynamic time 
scale, suggesting that the filling of the loop by the evapo- 
rated flare plasma may be the cause of the increase in the 
coronal column depth. It is straightforward to show that this 
is a distinct possibility (see the following chapter, Doschek 
etal, 1984). 

For the majority of observed bursts the most likely loca- 
tion of the microwave source is around the top of the loop, 
as discussed in Section 3.2.2. 1. However, the explanation 
for this observation is different than that presented above for 
hard X-rays. This is because the microwaves are the result 
of gyrosynchrotron emission by much higher energy elec- 
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trons (E — mc^) which, unlike the lower energy hard X-ray 
emitting electrons, cannot be confined easily. The confine- 
ment to the top of the loop of such high energy electrons 
by Coulomb collisions would require well above 10^^ 
cm"2. This is obviously unlikely and trapping in low den- 
sity (n < 10*® cm~^) loops or some other explanation is 
needed to explain the observed origin of the microwave emis- 
sion near the top of loops. 

Microwave-producing electrons are clearly present 
throughout the whole loop but the efficiency of gyrosynchro- 
tron emission varies throughout the loop and depends on the 
direction and strength of the magnetic field. The gyrosyn- 
chrotron radiation is strongest in the direction perpendicu- 
lar to the field line. This factor alone, i.e., assuming 
everything else is equal throughout the loop, will favor a 
maximum brightness located between the footpoints. For ex- 
ample, the simple loop of Figure 3.11a, when viewed from 
above, will be brightest at r = 0. Note that this aspect is, 
to the first order, independent of the viewing angle. 
However, the gyrosynchrotron emission depends also on the 
number of electrons and the strength of the magnetic field. 
At a given frequency v the emission from stronger field 
regions corresponds to lower harmonics of gyrofrequency 

The lower harmonics are produced by lower energy elec- 
trons, E - which are more numerous. Since the 

field strength only decreases with the height above the photo- 
sphere, emission from higher regions could dominate only 
if this variation is slow. Finally, the strength of the gyrosyn- 
chrotron emission is a function of the pitch angle distribu- 
tion of the electrons. Only if the pitch angle distribution of 
the accelerated electrons is broad will the emission from the 
coronal region be important. Otherwise, if electrons of kinetic 
energy E » mc^ are beamed along the field line, they re- 
main so throughout the corona and give little gyrosynchro- 
tron radiation. The pitch angle distribution becomes broader 
and the radiation efficiency increases only below the transi- 
tion region, i.e. , at the footpoints. Another factor which also 
aids the relative enhancement of emission from the coronal 
part of the loop is the fact that gyrosynchrotron radiation is 
subject to various absorption and suppression effects (cf. 
Ramaty and Petrosian 1972), all of which are more impor- 
tant in the chromosphere than in the corona. 

New observations of flare X-ray polarization by Tramiel 
et al. (1984), exhibit a low degree of linear polarization ( < 
5%). This result differs from earlier observation by Tindo 
and his collaborators (see e.g., Tindo et al., 1976). Such 
low polarization naturally will arise in the so-called thermal 
models or if the pitch angle distribution of the electrons is 
nearly isotropic. Earlier studies of nonthermal models (Haug 
1972, Langer and Petrosian 1977, Bai and Ramaty 1978), 
which ignored the scattering of the nonthermal electrons, con- 
sequently found a high degree of linear polarization for 
highly-b>eamed electrons. Even Brown’s (1972) calculation, 
which includes the scattering in an approximate fashion but 
does not include the curved loop geometry, has indicated a 


higher degree of polarization than observed. As a result it 
has been assumed that a low observed degree of polariza- 
tion is a strong evidence against the nonthermal models. 

Recent work by Leach, Emslie and Petrosian (1984) has 
examined the problem more fully, and does not confirm this 
expectation. The primary reason for a low degree of polar- 
ization is that no matter how strong the anisotropy of the 
injected particles, the average pitch angle distribution in the 
thick target is very broad (see Section 3. 2. 2.2). In addition, 
because the direction of the linear polarization varies along 
a loop, the integrated X-ray emission from the whole loop 
will have the characteristic of the emission from a highly 
broadened distribution and will show lower polarization than 
one would expect from the distribution of the injected elec- 
trons. Table 3.2 summarizes the effect of the spectral index 
5 and the pitch angle distribution dispersion parameter aj 
of the injected electrons, assuming the magnetic field to be 
uniform, and Figure 3.13 compares the observation of 
Tramiel et al. (1984) with three model calculations. 


Table 3.2 Maximum Percentage Polarization 
at 1 6 keV 
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A low degree of polarization can arise naturally from 
some nonthermal models. In particular, there is a strong ef- 
fect of spectral index on polarization; compare tabulated 
polarizations for various 6 values in Table 3.2. This is be- 
cause the degree of polarization increases rapidly as the pho- 
ton to electron energy ratio approaches unity (k/E — 1). For 
steep spectra (6 large) the bulk of photons with energy k are 
produced by electrons of only slightly larger energies 
(k/E ^ 1). However, for harder spectra (6 small) the con- 
tribution of higher energy electrons (i.e., those with k/E con- 
siderably less than unity) increases, causing a reduction in 
the polarization. 

Opposite polarization on both sides of the site of maxi- 
mum brightness (cf. Figure 3.3) is a natural consequence 
of the simple loop geometry we have been discussing here. 
The circular polarization is (e.g., Petrosian and McTiernan 
1983). 

Pc=cot 0 {[3(l+6)vin^/*']/[l+(Vi*^'^/^cosOT*/2}, ( 3 . 3 ) 
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Figure 3.13 A comparison between the polarization measured for flare 21 of 
Tramiel et at. (1984) and the polarization calculated for three of the models of 
Leach and Petrosian (1983). For the solid line aj = 0.4 and dfnB/di =0.0, the 
dashed line aj = 0.4 and dfnB/di = 1.5 x 10"®, and for the dotted line aj 
= 0.1 and d^nB/di = 1.5 x 10~®. A value of dfnB/dt = 1.5 x 10~® cor- 
responds to a twentyfold increase in the magnetic field strength from the top 
of the loop to the transition region at = 10’®cm^. For a flare at disk center 
the viewing angle would be 0°; for one on the solar limb it would be 90° (from 
Leach et aL, 1984). 


where 0 is the angle between the line of sight and the mag- 
netic field lines. As evident = 0 at S = 7t/2 and increases 
with opposite sense away from d = ttI2 where the total 
microwave intensity is at its maximum. The simple circular 
loop, therefore, agrees with the qualitative features of the 
observations. Alteration of the field geometry changes the 
variation of P^.. In particular, for a highly asymmetric loop 
(see e.g. , Figure 3.11b) it is possible that one sense of polar- 
ization becomes dominant, as often observed. Note that this 
is qualitatively similar to the field geometry proposed by 
Kundu and Vlahos (1979), though the details of their model 
and the one presented here are quite different. 

3.2.3 Combined Soft X-ray and Hard X-ray 
Spectra 

Analysis of the spatially integrated spectra from flares 
in the range 3 keV to 500 keV can yield much vital informa- 
tion regarding the separation of impulsive and thermal com- 
ponents of the flare. However, the interpretation is complex 
since several plasma regimes contribute to these spectra, and 
involve a range of atomic interaction processes. As the lower 


end, 3 keV to 7 keV, there are dominant line spectra of Fe 
XXVI, Fe XXV, Fe XXIV, Ca XIX, and in some cases Fe 
XXni to Fe XVn. Although these can be excited by ther- 
mal or nonthermal electrons, the existence of these species 
requires that the excitation occur within hot thermal plasma. 

Such hot thermal plasma also emits continuum radiation, 
which may extend up to photon energies of 10 keV, depend- 
ing on the temperature. Non-thermal electrons, or streams, 
passing through the hot thermal plasma deposit energy, and 
contribute an additional component of the continuum spec- 
trum. This process is normally represented by a thick target 
model. Power law electron distributions give rise to power 
law spectra. This deposition may take place within the hot 
plasma itself, or within cold dense plasma outside the hot 
region, as in the classic foot-point emission models. Although 
this alternative makes no detectable difference to the con- 
tinuum radiation emitted, it affects whether this mechanism 
contributes to the intensities of the above spectral lines. The 
range of mechanisms is thus large, and is further compli- 
cated since coverage of the region normally involves two or 
three spectrometers of quite different characteristics. 
However, it is clear that the ambiguity in interpretation is 
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greatly reduced by attempting a simultaneous analysis over 
a wide range of photon energies. 

3.2.3.I. Line Spectra 

The intensities of spectral lines, and their dependence on 
the plasma parameters, is determined by a range of atomic 
interaction quantities. For excitation of spectra from the 
ground state of the ion responsible, excitation and dielec- 
tronic capture rates are important. For helium and hydrogen- 
like ions, recent calculations by Bely-Dubau et al (1982a,b) 
and Dubau et al. (1981) are widely used. Further work is 
in progress on the lower ions Fe XXIII to Fe XVII. These 
spectra lie in the narrow range 1.85 to 1.91 A but cover 
a wide range of effective plasma temperatures. They are thus 
very valuable in extending the differential emission meas- 
ure analysis downwards below 10^ K, as well as in studying 
the transient ionization that occurs during rapid time varia- 
tions of flare plasma. Moreover, where the ground configu- 
ration of the ion is complex (Fe XIX to Fe XXII), the 
intensity ratio of some of the lines can be used as an elec- 
tron density diagnostic in the range 10‘^ to 10*^ cm‘^ This 
could be especially valuable in the early part of the impul- 
sive phase of hot flares. 

For many purposes ionization balance calculations are 
also required for interpreting the line spectra. Many sets are 
available in the literature. However, some new data are be- 
coming available from the solar flare spectra themselves. 
Antonucci et al. (1984, and further work in preparation) have 
used a study of the flare data from the BCS, Hinotori and 
P78-1 in order to derive ratios of Li-like to He-like ions in 
Ca and Fe as a function of temperature. The ionization 
balance calculations of Jacobs etal. {1911 ^ 1980) have been 
modified to take account of this work, and then used in some 
of the modeling described later. 

3.2.3.2 Transient Ionization 

It was earlier hoped that transient ionization effects would 
form an important diagnostic for the dynamic phases of 
flares. In this way, departures of the measured ionization 
ratios from their steady-state values would depend upon local 
electron density and the rate of change of temperature. It 
has since emerged that, due to the high electron densities, 
ionization is very close to steady-state over almost all phases 
of the flare for which statistically good data are obtainable. 
It has long been realized that the ratio of Li-like to He-like 
ions in flares is larger than theory would predict, by about 
a factor 2. This might be interpreted as a transient ionizing 
effect, were it not that it continues through all phases, in- 
cluding the flare decay. Koshelev and Kononov (1982) pro- 
posed a dynamic oscillating flare model in order to explain 
this effect by a continuous transient ionization condition. 
Doschek (1984) has attempted to examine such models on 
a more quantitative basis. He concluded that it is difficult 


to find an oscillating model which adequately predicts all the 
observed spectral intensities. A more likely explanation lies 
in a factor of 2 error in the ionization bdance models de- 
rived from theory. This assumption is implicit in the rederi- 
vation from observations by Antonucci et al (1984) 
mentioned above. 

The existence of real departures from steady-state ioni- 
zation are now being reported by some workers, who exa- 
mine critically the earliest impulsive phase data (see later 
in this section). If substantiated, this would provide a valua- 
ble measure of the density at this time. 

3.2.3.3 Hinotori Analysis 

The Hinotori team has attempted to analyze jointly data 
from two of the instruments; the flare monitor (FLM) and 
the hard X-ray monitor (HXM). The FLM covers the range 
5 to 20 keV with a resolution sufficient to separate line and 
continuum contributions. The HXM covers the range 18 keV 
to 400 keV. In an earlier analysis of FLM, Wantanabe (1984) 
found that while the continuum is best fitted by a power law 
in the early phase, it rapidly turns into a thermal spectrum. 
The time histories show that the thermal energy content is 
adequately represented by the time integral of the hard X-ray 
signal, indicating that the impulsive phase electrons may 
carry all the energy required for the flare. He concluded also 
that there is a strong tendency for all flares to produce a ther- 
mal component temperature of 2 x 10^ K and density 10'* 
cm~3. 

Tanaka et al (1984), in an attempt to understand the 
fluorescence origin of the Fe Ka line, have combined the 
spectra from FLM and HXM. The result. Figure 3. 14, shows 
a smooth blend between the power law behavior above 10 
keV with the thermal behavior below 10 keV. One is tempted 
to conclude that the plasma has a smooth distribution of elec- 
tron energies merging from streams to thermal around 10 
keV. This interpretation would be valid only if both compo- 
nents represented thin target emission. If the streams dissi- 
pate in thick target emission, then the transformation from 
photon spectrum to electron spectrum would be different in 
the two cases. 

3.2.3.4 SMM Analysis 

A method has been devised (Gabriel, Sherman, Bely- 
Dubau, Orwig and Schrijver, in preparation) to model the 
emitted spectrum based upon input in the form of a differ- 
ential emission measure distribution plus electron beams. The 
model computes the predicted fluxes from 8 spectral lines 
in the Bent Crystal Spectrometer (BCS), the 6 channels of 
the Hard X-ray Imaging Spectrometer (HXIS) and the 15 
channels of the Hard X-ray Burst Spectrometer (HXRBS). 
This then covers line emission from 10’ K to 10® K and con- 
tinuum energies from 3 keV to 260 keV, Beam energy dis- 
tributions can be arbitrary and need not follow power laws. 
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Figure 3.14 Combined spectra from the Hinotori FLM 
and HXM instruments, during the impulsive phase of 
a flare. 


Thin target emission by the beam within the hot plasma is 
evaluated as well as the normal thick target emission, but 
double counting of this radiation is avoided. An inversion 
technique is not used; it is rather a question of trying vari- 
ous distributions to see which best fit the data. In this way 
it is possible to take proper account of the varying accuracy 
of the data. For example, a satellite to resonance line ratio 
in one BCS channel should be good to * 10% , whereas ab- 
solute instrument calibration might be out by up to a factor 
of 2, although any such systematic errors must then always 
be present. 

The method has been applied to the flare of 1 822 TIT 29 
June 1980 at various times throughout the flare. The results, 
some of which are shown in Figure 3.15, are preliminary, 
but the following conclusions emerge: 

1 . A strong thermal component at 10^ K initially, rising 
to 2x10'' K later, is always present. 

2. This thermal component is not isothermal, but has a 
real width «0.5xl0^ K. 

3. At all times evaluated (even past the peak thermal 
emission) there is an additional high-energy compo- 
nent. If fitted with a thermal distribution, this is at 
T » 10* K and has an emission measure 3 or more 
orders of magnitude lower than the primary compo- 
nent. It can also be fitted adequately with a power- 
law electron beam. 

4. Early in the impulsive phase, the first thermal peak 
in the differential emission measure is smaller and at 


a lower temperature, while the higher energy com- 
ponent is relatively stronger. At the earliest time 
tested, it is very difficult to fit a thermal distribution 
to this feature. It is easier to fit a power law beam 
with a cut-off around 10 keV, with an electron spec- 
tral index of 5-7 on a thick target assumption. At this 
early stage, a good fit is obtained with the lines only 
if it is assumed that the ionization balance lags be- 
hind its steady-state value. 

Figure 3.15 (top) shows ratios of calculated to observed 
intensities for two times during the flare, early in the impul- 
sive phase and at the peak of the thermal phase. The 
anomalously high measured intensity of the Fe XXIV q line 
during the impulsive phase is an indication of the departure 
from ionization balance mentioned above. Figure 3.15 (bot- 
tom) shows the differential emission measures used in cal- 
culating these cases. 

It could be interpreted that the lower thermal peak 
represents the main flare loop being filled with evaporated 
chromospheric plasma, having a temperature gradient, and 
temperature stabilization produced by conduction and radi- 
ation cooling. The higher energy component is produced by 
either or both of thick-target bremsstrahlung from the 
primary fast electrons and a high temperature region iden- 
tifiable with the primary release. We thus have an indica- 
tion that for this flare, the energy input continues at a reduced 
level until well after the peak of the thermal phase. Such an 
interpretation would be complicated, but not invalidated, if 
the observed emission arises from two independent flaring 
regions. In any case, the high energy component, at its peak 
of intensity, is precisely identical to the normally-assumed 
signature of the impulsive phase, i.e. a signature of nonther- 
mal bremsstrahlung. 

3.2.4 Iron Ka Emission 

There has been great interest in the radiation mechanism 
of iron Ka, owing to its value as a diagnostic of nonthermal 
electron beams. The Ka lines at at 1.936 A and 1.940 A 
are radiated by inner-shell ls-2p transitions following the 
removal of a K-shell (Is) electron from iron ions Fel-FeXI. 
The inner-shell ionization occurs either by photoionization 
by X-rays of energy above the 7. 1 1 keV ionization threshold 
(fluorescence) or by electron impact ionization. Since the 
chromosphere is transparent for 7 keV photons, irradiation 
of cool photospheric iron by X-rays of E > 7 . 1 1 ke V could 
cause the emission (Neupert et aL , 1967). Electron impact 
ionization occurs most effectively for electrons with an 
energy of about 25 keV. It has been suggested that nonther- 
mal electrons associated with the hard X-ray burst excite Ka 
lines by collisionally ionizing the iron ions that are in lower 
ionization stages (Acton 1965, Phillips and Neupert 1973). 

Observations of Ka lines in the previous cycle, though 
performed with relatively low spectral and temporal resolu- 


3-17 



Ratio 



k r q w J q w 


FeXXV FeXXVI 


1 2 3 4 5 6 

I i I I I I 


1 2 3 4 5 6 7 8 

I I I I 1 I I i 


1822.40 

U.T. 


1825.00 

U.T. 




HXRBS 



PLUS POWER LAW \ 
ELECTRONS 

INDEX 5.3 ABOVE 25 kaV 
INDEX 7 BELOW 25 kaV 


PLUS POWER LAW ' 
ELECTRONS INDEX 7 


ALL THERMAL FIT 

^ ^ 


Figure 3.15 Top: Comparison of observed spectra and the model fit. Bottom: The differential 
emission measure used for the model fit shown above. 
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tion, have suggested fluorescence as a plausible mechanism 
(Doschek et al, 1971; Tomblin 1972). During the present 
maximum more reliable observations of iron Ka have been 
made with much increased spectral and temporal resolution 
and increased sensitivity. The P78-1 SOLFLEX (Feldman 
et al , 1980) and the SMM XRP (Culhane et al , 1981) ob- 
tained clearly resolved Kal and Ka2 lines using flat Ge and 
bent Ge crystals respectively. The Hinotori SOX (Tanaka 
et al , 1982) obtained barely resolved Ka lines and the K0 
line at 1.755A by a Si02 crystal, and the Tansei IV satellite 
(Tanaka 1980) obtained Ka by a LiF crystal which possessed 
the highest sensitivity. 

These observations have provided results which are in 
agreement. All flares show Ka emission which favors 
fluorescence mechanism almost throughout the flare. The 
fluorescence origin has been quantitatively established by 
Parmar et al . , (1984) from detailed comparison of the SMM 
observations of 40 large flares with Bai (1979)’s model of 
fluorescence. In Bai’s model, which used a Monte-Carlo 
technique to follow iron Ka photons resulting from the pho- 
tospheric absorption of X-rays, the Ka flux is evaluated as 
the products of three terms: integrated X-ray flux above 7. 1 1 
keV, fluorescent efficiency which depends on the X-ray 
source height and source temperature, and a term consider- 
ing the absorption and scattering of emergent Ka photons, 
which depends on the heliocentric angle of the flare. For each 
of the flares Parmar et al. , could show close similarity be- 
tween the Ka light curve and that of the E > 7.11 keV X- 
ray flux. In fact, the Ka light curves were quite different 
from the hard X-ray burst profile of E > 50 keV. A remark- 


able center-to-limb variation was found for the Ka to con- 
tinuum ratio for different flares at various heliocentric angles, 
and explained by Bai’s model (Figure 3.16). Referring to 
Bai’s model, the large dispersion in this ratio for flares at 
similar heliocentric angles is consistent with the variety of 
flare heights ranging from 0 to 0.05 solar radii for an iron 
abundance of 5.5 x 10~^ per H atom. Conversely, from 
the restriction for flare heights in short duration compact 
flares which were observed in the Skylab X-ray images 
(Pallavicini 1977), a plausible iron abundance has been esti- 
mated to be in the range from 5 x 10"^ to 6 x 10"^ per 
atom. 

While the Ka light curve shows gradual variations simi- 
lar to those of soft X-rays in most flares, transient Ka emis- 
sions in the very early phase of flare have also been reported 
in limited number of flares. Culhane et al. (1981) reported 
a Ka light curve showing impulsive behavior. During the 
15 seconds of the impulsive hard X-ray burst, the measured 
Ka flux was 2.5 ± 1.4 times that measured during the fol- 
lowing 15 seconds (Parmar et al . , 1984). However, this case 
has been considered rather exceptional in the SMM data. The 
LiF crystal on the Tansei IV satellite detected a sudden en- 
hancement of Ka emission that preceded an enhancement 
of the high temperature iron lines around 1.85 A in three 
flares (Tanaka 1980). This enhancement occurred in the ris- 
ing phase of the microwave burst at 17 GHz. The Hinotori 
observed similar, but less evident, cases in several (7-8) very 
impulsive flares. In these flares the Ka to Fe XXV resonance 
line ratio, which normally takes a fairly constant value of 
about 0.1, showed a very high value from 0.3 to 1.0 only 



Figure 3.16 The ratio of Ka flux to the E > 7. 11 keV continuum for 
various flares plotted against the heliocentric angle. The ratios pre- 
dicted by Bai's model for an iron abundance of 5.5 x 10"® per H 
atom, a temperature of 2.0 keV and three specific flare heights are 
shown as broken curves (Parmar et at., 1984). 
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in the initial phase, in close time coincidence with the hard 
X-ray burst. 

These examples may provide candidates for the electron 
impact Ka emission. However, Tanaka et al. (1984) have 
suggested that the observed excess emission is mainly ex- 
plained by fluorescence from the hard X-ray burst extend- 
ing in energy to 7 keV. The high resolution continuum 
spectra from the Hinotori have revealed that a single power- 
law distribution prevails from 70 keV to below 10 keV in 
the early phase of some impulsive bursts (see Figure 3,14). 
Because in the very impulsive flares the thermal X-ray flux 
originating from 20 million degree plasma at 7 keV is less 
than the flux of the extended power-law distribution, fluores- 
cence Ka shows impulsive behavior correlated with the hard 
X-ray burst. In the majority of flares, however, the thermal 
X-ray flux already dominates the 7 keV region when low- 
sensitivity crystal spectrometers can observe. One example 
of an impulsive event is shown in Figure 3.17. While the 
Ka light curve from the middle phase is explained by fluores- 
cence due to thermal X-rays, the Ka flux in the early phase 
much exceeds the predicted values. The excess fluxes are 
explained by fluorescence from the power-law X-ray fluxes 
assuming that the source is located at 0 km. The transition 
from power-law X-ray fluorescence occurs rather smoothly 
around the peak of the hard X-ray burst, as expected from 
the fact that the thermal X-ray flux increases most rapidly 
at that time, with its flux increasing approximately in propor- 
tion to the time integral of the hard X-ray flux (see, e.g., 
Tanaka et ai, 1982). However, the conclusion that all Ka 
emission is due to fluorescence depends on assigning a height 
of zero to the power-law source. If this source is located high 
up in the corona, as has been observed in some gradual bursts 
observed by Hinotori, then an excess Ka flux remains to be 
explained, presumably by electron impact. In the impulsive 
bursts like this case, however, this seems unlikely. 

This result may be compared with the thick target theory 
of the electron impact Ka developed by Emslie, Phillips, and 
Dennis (1985). The electron impact Ka yields predicted by 
this theory are always smaller than the Ka yields for fluores- 
cence due to the power-law X-rays extending down to 7 keV, 
provided that the spectral index is larger than 3.5. There- 
fore, in most cases the electron impact Ka emission is veiled 
under the fluorescence Ka, even if the electron beam does 
exist. In conclusion, the electron impact Ka appears to be 
not necessary to explain the observations. On the other hand 
an evidence of low energy (below 10 keV) power-law pho- 
tons may raise another problem related to the flare energetics 
because the energy included in the low energy electrons as 
derived from the thick target hard X-ray model is overwhelm- 
ingly large, probably exceeding the energy radiated in the 
UV and X-rays regions (Tanaka et al, 1984). 

3.2.5 Ultraviolet and Hard X-ray Emission 

The spatial, temporal and energetic relationships between 
impulsive ultraviolet emission and hard X-ray emission have 
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Figure 3.17 Ka light curve for an impulsive burst of 
1 981 August 1 1 (Observed by Hinotori), together with 
the Ka intensities predicted from fluorescence due to 
the thermal X-ray flux as derived from the iron line 
spectra (thick line) and those predicted from fluores- 
cence due to the observed X-ray continuum which 
showed a power-law spectrum down to 7 keV 
(dashed line). Bars indicate ± a range of the observed 
counts. Time profiles of the hard X-ray burst, Fe XXV 
resonance line intensity and electron temperature de- 
rived from the line spectrum are also shown. 

attracted considerable interest as tests of flare particle trans- 
port models. Kane et ai (1980) reviewed the status of this 
field in the Skylab era. Two quite different UV observational 
approaches have been taken. Below, in Section 3.2,5. 1, we 
review and interpret observations in the wavelength band 
from 10 to 1030 A. This band contains emission lines and 
continua that originate at a wide range of temperatures, from 
lO^-lO"^ K. We then turn to recent SMM observations, which 
are more straightforward to interpret— observations of the 
transition-zone line 0 V 1371 A , discussed in Section 3.2. 5. 2, 
and their intepretation in terms of both thermal and nonther- 
mal energy-transport models. 

3.2.5.1 10 - 1030A and Hard X-ray Emission 

An extensive set of calibrated concurrent measurements 
of impulsive hard X-ray and integrated 10-1030 A (hence- 
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forth EUV) emission has been obtained by Kane and 
Donnelly (1971), Donnelly and Kane (1978) and Kane, Frost 
and Donnelly (1979). The EUV observations have no spec- 
tral or spatial resolution, but good temporal resolution 
( ~ 1 s). Although these measurements were made indirectly 
using the response of the ionosphere, subsequent direct 
meaurements from spacecraft by Horan, Kreplin and Fritz 
(1982) have fully validated the indirect method. The hard 
X-ray observations have no spatial resolution, but have suffi- 
cient spectral resolution to allow a power-law spectral fit, 
and also have good temporal resolution. All studies of the 
correlation between the EUV and hard X-ray bursts have 
found coincidence of peaks in the EUV flux and X-ray flux 
to within the time resolution of the instrument. This observa- 
tion rules out certain types of thermal models (see Section 
3.2.5.2b). 

The observational datum of primary interest for each im- 
pulsive burst is the EUV/HXR ratio at the peak, where HXR 
is the flux in hard X-rays ( > 10 keV), both measured at the 
top of the earth’s atmosphere. Data for many impulsive bursts 
are shown in Figure 3.18. Donnelly and Kane (1978) and 
Kane, Frost, and Donnelly (1979) concluded that the peak 
EUV/HXR ratio is constant, rejecting nonlinearity in the peak 
EUV/HXR relationship seen in large flares as a selection ef- 
fect. On the other hand, McClymont, Canfield and Fisher 
(1984) concluded that the peak EUV/HXR ratio is not con- 
stant, but rather the peak EUV a HXR^^^ jssyg arises 
because the largest flares in Figure 3. 18, for which the peak 
EUV/HXR ratio is smallest, tend to have been observed by 
a single spacecraft, TDl-A. 

Until recently there has been general agreement that the 
observed EUV/HXR ratio is about 30 times smaller than the 
value predicted by typical electron beam fluxes in the thick- 
target nonthermal model, in which all the X-ray emission 
is produced by nonthermal bremsstrahlung and all the EUV 
emission by thermal emission from the collisionally-heated 
transition region and chromosphere (Donnelly and Kane 
(1978), Emslie, Brown and Donnelly (1978), Kane, Frost, 
and Donnelly (1979)). Donnelly and Kane (1978) suggested 
that the discrepancy could be explained by the “partial 
precipitation” model, in which the bulk of the X-ray emis- 
sion is from electrons trapped in the corona and only a small 
fraction of the nonthermal electrons reach the chromosphere. 
However, Emslie, Brown and Donnelly (1978) pointed out 
that trapping of nonthermal electrons in the corona could not 
cut down the EUV flux sufficiently because electrons are 
scattered out of the trap onto the chromosphere after losing 
only roughly half of their energy (Melrose and Brown, 1976). 
Emslie, Brown, and Donnelly (1978) suggested that the ther- 
mal hard X-ray model (e.g.. Brown 1973a) may be more 
appropriate. 

The peak EUV/HXR ratio also shows a strong center- 
to-limb dependence. Donnelly and Kane (1978) showed that 
this dependence could be explained if the EUV emission 
comes from a source embedded in the chromosphere. They 



Figure 3.18 Observed retationship between the peak 
10-1030 A EUV/HXR ratio and the HXR energy flux 
above 10 keV measured at 1 AU, at peaks of impul- 
sive bursts, from Kane and Donnelly (1971), Donnelly 
and Kane (1978), and Kane, Frost, and Donnelly 
(1979), McClymont, Canfield, and Fisher (1984) have 
multiplied the observations in the 20-70® central 
meridian distance (CMD) range by 2.5, as an approx- 
imate correction for the effect of limb darkening. 
Lines: Predictions of the thick-target nonthermal 
model, for two different values of the photon spec- 
tral index 7, and a universal beam cross-section of 3 
X 10’® cm^. From McClymont, Canfield, and Fisher 
(1984). 


explained the dependence satisfactorily by a model in which 
the EUV emission comes from the bottoms of cylindrical 
wells of two different sizes. The idea that the EUV flare 
source is below the normal nonflaring chromosphere is com- 
patible with observations of Ha spectra discussed in §n.G 
below. 

McClymont, Canfield and Fisher (1984) have looked at 
another aspect of the predictions of the thick-target model. 
Following previous authors (Brown 1973b, Lin and Hudson 
1976, Emslie, Brown and Donnelly 1978), they calculated 
the depth of the flare transition region as a function of the 
electron beam flux per unit area. Considering only coUisional 
energy losses, as in previous work, they find that the increase 
in depth of the transition region with increasing beam flux 
causes the EUV/HXR ratio to fall as the flux increases, 

EUV/HXR oc (HXR/Area)-<^-2)/(«-^2), (3.4) 
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where the electron spectral index 5 is related to the hard X-ray 
spectral index 7, 5 = 7 + 1. The physical reason for the 
ratio to decrease with increasing flux is that the mass frac- 
tion of the beam-heated column that is filled with chro- 
mospheric material, which emits much of the EUV, decreases 
as flux increases. For the most commonly observed spectral 
indices, 7 « 4, the predicted slope of the peak EUV/HXR 
ratio as a function of hard X-ray flux is in excellent agree- 
ment with the slope that they find in the observational data, 
if it is assumed that the beam cross-section is constant in 
all flares y i.e. only the electron energy input per unit area 
varies. Second, as shown in Figure 3.18, both the slope and 
the absolute value of the peak EUV/HXR ratio vs. HXR 
curve are correctly predicted if one assumes that the beam 
cross-section is of order 3x10^^ cm^, comparable to typ- 
ically observed chromospheric flare kernel areas. 

Were it not for the fact that extremely high beam fluxes 
are implied, one would conclude that the thick-target elec- 
tron heating model, with a universal beam cross-section, ex- 
plains the observations very well. With a 20 keV cutoff 
energy, the beam fluxes range from about 10*® ergs cm"^ 
s"* for the smallest flares to 10*^ ergs cm"^ s"* for the larg- 
est. The majority of the observed flares have beam fluxes 
greater than 10** ergs cm'^ s"*, the commonly accepted 
upper limit for return current stability. 

How, then, can the peak EUV/HXR ratio be understood, 
if beam fluxes greater than about 10** erg cm"^ s~* are dis- 
allowed? McClymont, Canfield and Fisher (1984) considered 
three processes which decrease the EUV/HXR ratio predicted 
by the thick target model. Firstly, if the preflare loop, in 
which electrons are acclerated, is hotter (and therefore 
denser) than a typical preflare loop—a state that might be 
inferred from the frequent observation of a soft X-ray precur- 
sor to the impulsive phase— the bulk of the electron beam 
energy may be deposited in the corona. They found that the 
drop in peak EUV/HXR ratio could be explained by this 
hypothesis, although rather high preflare temperatures 
( > 10^ K), and correspondingly high densities, were neces- 
sary to obtain the very low peak EUV/HXR ratio observed 
in large flares. Secondly, thermal hard X-rays may be im- 
portant, as investigated previously. Thirdly, the effects of 
reverse (return) currents have not been considered fully in 
calculations of this type. Due to reverse currents, concen- 
tration of heating at the foot of the transition region may 
greatly increase the ability of electron beam heating to raise 
chromospheric plasma to coronal temperature, compared to 
pure collisional heating, thus reducing the emission meas- 
ure in the temperature range where the EUV emission 
originates. No quantitative studies of this effect have yet been 
done. 

3.2.S.2 OV and Hard X-ray Emission 

a. Observations. Observations of individual EUV lines 
are potentially more powerful than measurements of the total 


10-1030 A EUV emission, discussed above. One such line, 
well observed by SMM, is the OV 1371 A line. The SMM 
observations have demonstrated a close observational rela- 
tionship between this line, formed in the transition region 
at a temperature of about 250,000 K, and hard X-rays. 
Progress has been made at defining temporal, spatial, ener- 
getic, and morphological relationships. 

The temporal relationship between HXR and OV was 
studied in detail by Woodgate et al (1983). In this work they 
studied the three SMM flares for which UV data were avail- 
able on a time scale of 1.3 s, and HXR data (from 25—100 
keV) were available every 0.128 s. Each flare displayed at 
least five impulsive bursts in both HXR and OV that were 
simultaneous to within the 1.3 s time resolution. Measure- 
ment of time lags for the peaks (tQy - tj^) yielded an aver- 
age of 0.3 ± 0.5 s. Thus, their results showed that peaks 
in emission that occurred in both OV and HXR were simul- 
taneous to well within the time resolution. 

The Woodgate et al (1983) results also showed that there 
was a very good occurrence correlation between OV and 
HXR peaks. For the 16 HXR peaks there were 14 OV peaks 
(80%) and for 19 OV peaks there were 14 HXR peaks (74%). 
The correspondence is obviously quite good “with the 
primary difference being due to some small UV peaks late 
in the impulsive phase that do not appear in the hard X-rays’ ’ . 
These small OV peaks may be associated with HXR bursts 
that are below the threshold of detectability, but may also 
indicate additional sources of heating. 

Poland et al, (1984) have studied the relation between 
energy emitted in HXR vs. OV as a function of time for 
several flares. The results for the flare of 2 November 1980, 
whose light curves are shown in Figure 3.19 are pre- 
sented in Figure 3.20, In the latter figure the boxes represent 
data points before flare maximum, the diamonds are within 
5 s of flare maximum and the plusses are after maximum. 
It can be seen that there is a well-defined relation through- 
out the flare between the HXR and OV emission. Also, the 
premaximum and postmaximum emissions lie on almost the 
same line, with the postmaximum OV remaining slightly 
brighter than the premaximum emission for the same HXR 
emission. Ten flares were studied in this manner, with all 
showing similar results. Most of the flares show the rise 
phase and fall phase occurring on approximately the same 
line. However, some do show a significant deviation from 
this relationship in that OV emission is brighter during the 
decay phase than during the rise phase. There is also a sig- 
nificant variation in OV emission from flare to flare, but, 
this may in part be due to the variation in total background 
in OV from one flare to the next. When these data are 
plotted with the preflare background removed there is still 
a large variation in the OV intensity from one flare to the 
next. 

The most important point to be made from these obser- 
vations is that HXR and OV emission rise and fall together, 
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Figure 3.19 The upper curve shows emission in OV for the entire flare 
as a function of time (actual values are 10*® times larger than shown). 
The lower curve shows emission in HXR above 25 keV for the entire 
flare as a function of time (actual values are 3.84 x 10**’ times larger 
than shown). Noise errors for a single measurement are shown at two 
different times for each plot. 



Figure 3.20 HXR emission as a function of OV emission in ergs s"’ for the 
flare shown in Figure 3. 1 9. Squares represent data obtained before flare maxi- 
mum, diamonds are within 5 s of maximum, and plusses are after flare maxi- 
mum. Only every third point has been plotted for clarity of the figure. 
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with the rise and fall phases usually on approximately the 
same line, independent of time, during the entire event. How 
can this be understood? It suggests that the physical condi- 
tions in the transition region are directly related to the produc- 
tion of HXR, and independent of those flare quantities that 
are known to vary on gradual phase timescales, such as the 
temperature or emission measure of the thermal X-ray 
plasma. It is straightforward to show that on gradual phase 
timescales ( ^ 1 m) the emission measure at OV tempera- 
tures in a conductively dominated transition region is propor- 
tional to the transition region pressure, so the observation 
of relatively little intensity variation during a flare such as 
the one shown in Figure 3.20 implies that the transition region 
pressure in the emitting regions varies very little in the course 
of the observations (a period — 1 m). If the model of an elec- 
tron beam in a single flare loop is to explain both the impul- 
sive and gradual emissions, then apparently pressure 
equilibrium is either never reached or is substantially un- 
altered by the impulsive event. On the other hand, it may 
be that the impulsive phenomena arise in a different volume 
than the gradual phenomena; we explore the arguments on 
this issue below, pro and con. 

To within SMM’s ability to spatially resolve emitting 
regions in HXR and UV, the two emissions arise from the 
same location on the sun. However, one must bear in mind 
that the UV observations are obtained at only 3 " x 3 " spatial 
resolution, at best, and more usually 10"x 10". X-rays be- 
tween 3.5 and 30 keV are observed with a resolution of 
7 "x7 " (van Beek et al, , 1980), at best. Harder X-rays, above 
30 keV, are observed without spatial resolution. Hence, our 
knowledge of the cospatiality of these emissions is strongly 
constrained by observational limitations. 

Spatially resolved observational information on the 
cospatiality of UV and HXR emission is limited to just a few 
flares. In a study of the June 29, 1980 flare Poland et al, 
(1982) showed that positions and changes in positions of UV 
transition region emission corresponded to positions and 
changes in position in 3.5 to 30 keV HXIS emission. These 
correspondences, together with the temporal correspondences 
between OV and higher-energy HXRBS hard X-ray light 
curves strongly suggest their cospatiality in this flare. Duijve- 
man et al, (1982), in a study of the November 5, 1980 flare, 
reported cospatial HXR and UV bursts in one of the flare 
footpoints. Machado etal, (1982) found similar cospatiality 
in the April 10, 1980 flare. Although these two latter papers 
did not study the UV/HXR correlation in detail, the HXIS 
data quality is better than for the June 29 event. 

In another SMM flare observation, Cheng et al, (1981) 
examined the SilV and OIV UV emission with a resolution 
of 4"x4" and found that each individual burst was seen in 
both UV and spatially unresolved hard X-rays, and that each 
component U V spike originated from a separate discrete flar- 
ing point. The density at these points was determined to be 
in the range of 5x10*^ to 10^^ cm'^. There is an uncer- 
tainty in using this ratio since the two lines are formed at 


significantly different temperatures; the reliability issue is 
discussed in Cheng et al, (1982). 

The close time-independent functional relationship be- 
tween OV and HXR, shown in Figure 3.20, motivates us 
to ask what is the obervational evidence that indicates the 
relationship between the volumes in which impulsive phase 
emissions and gradual phase emissions arise? Two argu- 
ments suggest that impulsive phase emissions and gradual 
phase emissions arise in separate volumes. Firstly, it is 
known that in X-rays the emission is impulsive only for ener- 
gies above approximately 5 keV; for lower X-ray energies 
the emission does not vary impulsively. For transition region 
emission, flare enhancements become increasingly impulsive 
above 10"* K. This impulsiveness reaches a maximum near 
2 X 10^ K, decreases rapidly at 3 x 10^ K, but remains sig- 
nificant to approximately 2 x 10^ K, above which the emis- 
sion is no longer impulsive (Donnelly and Hall 1973). 
Secondly, Widing and Hiei (1984) showed that impulsive 
emission in one flare arose from a single loop, while grad- 
ual phase emission arose from a separate location. The im- 
pulsive phase in this flare was seen in hard X-rays (20-30 
keV) and in transition zone lines between approximately 
80,000 K and 2 X 10^ K. The gradual phase emission from 
this flare, seen in soft X-rays and UV lines above 2x10^ 
K, was observed to arise from a separate region or loop near 
the impulsive loop. 

On the other hand, several arguments can be made against 
the assertion that impulsive and gradual emissions occur in 
separate volumes. Duijveman, Somov, and Spektor (1983) 
showed that one can see impulsive variations in bands as low 
as the 3.5-8 keV band of HXIS. Also, there is a well-known 
correlation between the amount of hard X-ray emission and 
the slope of the soft X-ray light curve in a typical flare soft 
X-ray emitting line (see, e.g. , Tanaka et al, , 1982a, Machado 
(1983). Finally, in some cases impulsive-phase phenomena 
occur in loops that already produce considerable thermal 
emission (see, e.g., Duijveman, Somov, and Spektor 1983); 
it is quite possible that the impulsive soft X-ray component 
is masked within the overall emission, as a spectrum such 
as Figure 3.14 would suggest. Only with substantially bet- 
ter spatial resolution will this issue be decided in a totally 
satisfactory manner. 

In summary, the impulsive phase observations clearly 
show that there is a strong physic^ link between the processes 
that form HXR and transition zone and chromospheric emis- 
sion. This link was established for flares in general by the 
relation in energies of peak emission found by Kane and 
Donnelly (1971), see Section 3.2.5. 1. The simultaneity of 
individual peaks in HXR and OV in several flares found by 
Woodgate et al, (1983), and the relation in energy emitted 
during several flares as described by Poland et al (1984), 
places observational limit on the possible links. There can- 
not be a significant time delay (i.e. greater than 1 s) between 
the different excitations, and the energies exciting the emis- 
sions are functionally related. That the emissions occur from 
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the same location to within 10 " supports the concept of their 
being physically related, but this not a serious constraint since 
this large distance could allow quite separate locations in the 
same loop or even different loops. 

b. Interpretation. Emslie and Nagai (1984) have studied 
the response of a coronal loop to both heating and nonther- 
mal electron acceleration at the loop top, in order to deter- 
mine the resulting structure of the transition region plasma 
(T » 10^±* K), and so the emission measure in, and inten- 
sity of, optically thin emissions formed in this region. The 
key quantity in the modeling of Emslie and Nagai (1984) is 
the logarithmic differential emission measure 

KT) = n2 dz/dfnT (3.5) 

where n is plasma density, z vertical distance, and T elec- 
tron temperature. Under an optically thin approximation for 
the radiative losses from the gas, this quantity is directly 
proportional to the intensity in a given line, the other fac- 
tors being atomic in nature and therefore independent of the 
structure of the atmosphere (e.g., Pottasch 1964). Thus de- 
termination of the evolution of {(T) with time t gives a model 
prediction of the shape of the intensity versus time profile 
of the line emission being considered. 

Figures 3.21 and 3.22 show the form of $(T) for various 
times since the onset of energy input into the flare loop for 
the conductively heated and electron-heated models, respec- 
tively. It is readily apparent from these figures that the evo- 
lution of J is markedly different in the two heating scenarios. 
Let us consider first the conductively-heated model, in which 
the rate of energy deposition varies linearly with time, up 
to a maximum of 2 x 10^® erg cm"^ s"‘ after 30 seconds, 
and in which the energy deposition is spatially distributed 
over a Gaussian profile, symmetric about the loop apex, with 
dispersion a = 2 x 10* cm. In this model, there is virtually 
no change in $ for the first few seconds, reflecting the fact 
that the thermal energy has not yet had time to propagate 
to the transition region layers. As soon as the heat front ar- 
rives (at t = 7 s for the parameters considered), a decrease 
in f results, due to the strong steepening of the temperature 
gradients produced by the large conductive flux F^( a T^^^ 
dT/dz), which evidently overcomes the increase in the n^ fac- 
tor (Equation (5)) in the determination of f . A few seconds 
later, however, the transition region plasma overheats to 
coronal temperatures, and the gas at the relevant tempera- 
tures ( — 10^ K) is now situated much deeper in the at- 
mosphere, with a corresponding larger density. The density 
scale height is sufficiently small in the preflare chromospheric 
layers that the conduction front, which moves at speeds of 
order of the local ion sound speed c^ = (kT/mp)^^^^ where 
k is Boltzmann’s constant and nip the proton mass, covers 
the distance to the “new” transition region in a very short 
time. Thus $ abruptly rises by several orders of magnitude 
in a second or two. This predicted trend in $(10^ K) (and 



Figure 3.21 The logarithmic differential emission 
measure (DEM) equation (3.5), for various times (in 
seconds) in the impulsive phase, for the conductively 
heated model. The dashed line shows the initial (t = 0) 
profile. Note the behavior at transition region 
temperatures—constant for the first few seconds, fol- 
lowed by a sudden decrease and an equally abrupt in- 
crease. The peak at around T = 1 0^ K on the 8 s curve 
IS rcalj it ccrrespcnds to the density enhancement that 
causes the large increase in E, at transition region tem- 
peratures in the 10 s curve. 


so in the intensity of optically thin EUV emissions) is simply 
not reflected in the observations above, implying that such 
an energy mechanism is not a viable possibility. 

On the other hand, for the thick target electron-heated 
model, in which energy is distributed throughout the flare 
loop with an identical time structure but with the broader 
vertical spatial structure resulting from consideration of the 
collisional dynamics of an electron beam in a relatively cool 
target (see Emslie 1978), the strong steepening of tempera- 
ture gradients associated with a relatively local energy 
released region does not happen. Instead, the dominant ef- 
fect controlling the behavior of $ is the progressive over- 
heating of denser and denser layers to coronal temperatures. 
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Figure 3.22 The logarithmic differential emission 
measure (DEM) equation (3. 5), for various times in 
the impulsive phase, for the parallel-injection electron- 
heated model. Note the steady Increase around 
transition-zone temperatures (T — 10® K). 


resulting in a monotonic rise of the density at T — 10^ K 
and so in ^ (10^ K). 

This correlation of EUV intensity with energy input [as 
measured by the instantaneous hard X-ray flux, which is 
roughly proportional to injected electron energy (Brown 
1971)] is in accord with the above observations (see, e.g., 
Poland et aL , 1982, 1984; Woodgate et al . , 1983). Indeed, 
Woodgate et al (1983) have used the close synchronism of 
EUV and hard X-ray emissions to rule out conductively 
heated models such as that discussed above. The results of 
Emsiie and Nagai further support this rejection. However, 
they also demonstrate that the observed behavior is predicted 
by the thick target electron heated model. On reflection, it 
is by no means obvious that this should be so. The observed 
increase in intensity is only an order of magnitude or so even 
in large events (Poland et al, 1984), a factor well below 
that predicted by the n^ factor in Equation (5), since n rises 
by over two orders of magnitude as a result of flare heating 
of the overlying layers (e.g. , Machado et al . , 1980, Emsiie 
and Nagai 1984). Thus the increase in the temperature gradi- 
ent factor must also play a role, and it is somewhat remark- 
able that these two large effects in opposite directions 
combine to produce a relatively small overall increase in 


as is observed. The underlying physical reasons for this be- 
havior are not fully clear at present, and merit further in- 
vestigation. 

3.2.6 White Light Emission 

We define white light flares as those detected in the opti- 
cal continuum; thus their number is limited by observational 
considerations. In the last few years patrols at the National 
Solar Observatory, Sacramento Peak (NSO/SP) and Big Bear 
Solar Observatory (BBSO) have greatly increased the known 
number [see catalog of white light flares by Neidig and Oliver 
(1983)]. It is probable that all major flares have detectable 
continuum in the blue, although the number detected at longer 
wavelengths is still small. Thus the term “white light flares” 
simply means highly energetic flares. Our general motiva- 
tion in studying white light flares is based on the great energy 
involved, the problem in transporting energy to the low at- 
mospheric heights of significant continuum opacity, and a 
general fascination with intense flares. An important spe- 
cial interest here is in the nature of the continuum emission 
and how it may be produced. 

3.2.6.1 Morphology 

Our knowledge of the morphology of white light flares 
is limited by the small number of events for which good data 
are available. Nevertheless, four kinds of white light flare 
have been recorded: (1) Bright footpoints, usually at the edge 
of the sunspot umbra (Zirin and Neidig 1981); (2) Short- 
lived flashes (Zirin and Tanaka 1981) associated with im- 
pulsive X-ray spikes, up to a few arcsec diameter and 20 
s lifetime, usually located along the neutral line of a delta 
configuration; (3) Moving fronts (Machado and Rust 1974; 
Zirin and Tanaka 1973) occurring after the impulsive phase 
and possibly associated with the hot thermal plasma; (4) Fixed 
bright points or ribbons (Zirin and Neidig 1981 ; Neidig and 
Cliver 1983), possibly associated with the hot thermal phase. 

Except for the last, all the manifestations are moving, 
sometimes very fast, or jumping from point to point. The 
white light flare footpoints commonly appear at a number 
of locations and times throughout the whole event; presum- 
ably they are associated with the main bundle of field lines 
on which flare energy input occurs. Almost always at least 
one footpoint occurs in the penumbra, usually at the edge 
of the umbra, but other bright points may occur outside the 
sunspot, as do white light flare ribbons. The area of white 
light flare emission averages — 6x lO^"^ cm^ at flare maxi- 
mum, with a duration ordinarily less than 10 minutes (Neidig 
and Cliver 1983). 

The white-light kernels and fronts roughly match the 
emission observed in He D3 and the wings of Ha. The He 
D3 emission, which may cover a somewhat larger area than 
the associated white-light emission, becomes especially bright 
in the post-impulsive phase (for example, the July 1 flare 


3-26 





(Zirin and Neidig 1981)). Probably all white light flares are 
accompanied by D3 emission. 

Very little information is available regarding the white 
light flare relationship with spatially-resolved, high-energy 
emissions. In the June 6, 1982 flare (Ohki et al , in prepa- 
ration) the Hinotori hard X-ray and soft X-ray images gener- 
ally show spatial and temporal agreement with kernels of 
emission at X3862 A . Videotapes show small kernels a few 
arcsec across, with changes occurring on time scales of a 
few seconds. 

3.2.6.2 The Optical Spectrum 

a. Observations, Spectral data on white light flares are 
usually obtained either from a series of broad band filter- 
gams at several widely spaced points in the spectrum (e.g., 
Zirin and Neidig 1981) or from spectrograms with disper- 
sions of several A /mm (e.g., Neidig 1983). In filtergrams 
the variation of flare intensity with wavelength is affected 
by changes in seeing between the successive images, or even 
by temporal changes in the flare itself. Fast switching with 
video measurements removes some of these problems, and 
the seeing may be evaluated. 

Another problem in interpreting the filtergram measure- 
ments is the possible effect of lines in the blue. Although 
the Balmer line emission can be minimized by using, for ex- 
ample, a narrow (15 A) filter at 3862 A (midway between 
Hg and H9), there still remains the possibility of emission 
by numerous photospheric lines (Hiei 1982). 

The spectrograms, although able to record a large ex- 
panse of the spectrum in a single exposure, are typically ob- 
tained with lower spatial resolution, with the result that the 
flare emission may be highly diluted or the most intense ker- 
nels may not fall on the slit. Judging by the intensities ob- 
served, the spectra obtained so far have not recorded the most 
intense kernels. Intensities up to twice the continuum at visi- 
ble wavelengths have been recorded with filters, while the 
spectral data have only revealed enhancements of 10-20% 
above the continuum. In addition, the variation of intensity 
with wavelength, as determined from spectrograms, may 
suffer from systematic errors due to chromatic aberration 
of atmospheric or instrumental origin. Thus while the gross 
features of the white light flare spectrum, the line intensi- 
ties, and the Balmer jump can reliably be obtained from spec- 
tra, the slope of the continuum thus obtained, which might 
be used to derive the flare temperature or to identify the emis- 
sion mechanism, should be treated with caution. 

The filter measurements, as well as some of the spectro- 
grams, show a fairly flat continuum in the visible with an 
increase below 4000 A . Superposed on this there is often 
a relatively strong Balmer continuum. Figure 3.23 combines 
filter data on several flares; the data for July 1 and April 
24 are from photographic measurements; that for May 9 is 
from high speed video measurements. Figure 3.24 shows 
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Figure 3,23 Broad-band filtergram measurements of 
the intensity enhancement [{flare — background)/ 
background] at flare maximum for several white light 
flares. Measurements at 3610, 4275, 4957, 5650 
and 6203 A are from NSO/SP; 3400, 3862, 4642, 
and 5400 A are from BBSO. 

spectrographic data in the vicinity of the Balmer jump for 
three representative flares. 

The Balmer continuum has been identified in spectro- 
grams of six events (Hiei 1982; Hiei et aL, 1982; Neidig 
and Wiborg 1984; Donati-Falchi et al, 1984), although 
several cases are known (e.g., Machado and Rust 1984; 
Boyer et al . , 1985) in which a short wavelength continuum 
was present, but without a measurable Balmer jump. A 
Paschen jump was detected in the spectrum of the 24 April 
white light flare (Neidig and Wiborg 1984); unfortunately, 
the latter flare is the only known event with a infrared con- 
tinuum sufficiently strong to permit a conclusive search for 
the Paschen jump. 

In addition to the large intensities observed in the Balmer 
continuum, the increase in the flare contrast in the blue, es- 
pecially below 4000 A , can be quite striking in some cases 
(Zirin 1980, Zirin and Neidig 1981). This “blue continu- 
um” is probably related to the bluish color noted in a num- 
ber of visual observations of white light flares, although the 
cause for it is not understood. An increase in the flare con- 
trast could be expected simply as a result in the decline in 
intensity of the solar background at short wavelengths (this 
would be particularly evident in the broad band filtergrams 
where the solar background is blanketed by lines). On the 
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Figure 3.24 Spectrographic measurements for three white light flares: Apr. 
24, 1981 (Neidig 1983); Sept. 10, 1974 (Hiei 1982); Aug. 7, 1972 
(Machado and Rust 1 974). lo refers to the intensity in the bright windows 
of the spectrum adjacent to the flare. 


other hand, intense brightenings (AI/Iq 1) have been ob- 
served in filtergrams at 3835 or 3862 A when the cor- 
responding emission at longer wavelengths was weak or 
absent (Zirin 1980). The later observations indicate differ- 
ences in contrast too large to be attributed entirely to differ- 
ences in the intensity of the solar background. The question 
as to whether the blue continuum is truly anomalous, aris- 
ing, perhaps, from some unidentified source of opacity or 
from photospheric lines, cannot be answered until more pre- 
cise observations of the spectrum become available. 

b. Interpretation.The principal task in interpreting the 
spectral observations is to identify the emission mechanism; 
this, in turn, determines the general atmospheric regime of 
the white light flare phenomena. The spectra in Figure 3.23 
suggest a strong variation with wavelength in the opacity of 
the white light flare source, indicating an origin in layers 
that are not optically deep. The Balmer continuum (Figure 
3.24), as well as the bright He D3 emission associated with 
white light flares, is characteristic of chromospheric densi- 
ties (electron density 1-5x10*^ cm"^; see Machado and 
Rust 1974; Neidig and Wiborg 1984; Zirin 1983). Neverthe- 
less, it is by no means certain that the emission underlying 
the Balmer continuum and extending to longer wavelengths 
is due to emission in general (the 24 April flare, which 
shows a Paschen jump, is a good candidate for an event domi- 
nated in the visible and infrared by emission). Hiei 
(1982) has pointed out that the continuum extending to longer 
wavelengths might be attributable to H" emission, possibly 
originating in the photosphere; this hypothesis remains un- 


tested because information on the Paschen jump is not avail- 
able for the flares studied by Hiei. 

White light flares in which the Balmer jump is absent 
(e.g., August 7, in Figure 3.24) might suggest either an op- 
tically deep source with small temperature enhancement or 
an optically thin, very hot source with opacity dominated 
by free-free absorption. The latter possibility, however, 
might require a prohibitively large emission measure at high 
temperature (n^^ Az > 10^^ cm-^ at 10^ K, in order to 
produce a 10% enhancement at 4QOO A). The spectrum of 
the August 7 flare was studied by Machado and Rust (1974) 
who concluded that the emission was Paschen continuum 
originating in the flare chromosphere. The latter interpreta- 
tion, of course, would require Balmer continuum to be 
present in this flare, and therefore, unless the absence of the 
Balmer jump can be explained, a deeper, low-temperature 
source with H" opacity might be indicated. In any case the 
height of the source in the August 7 flare would have to be 
greater than 300 km above T5000 = 1, as Machado and Rust 
found a probable absence of emission in lines formed below 
that level. With regard to a deep photospheric origin for white 
light flares in general, it is worth noting that strong rever- 
sals in the cores of photospheric absorption lines have not 
been observed in white light flare spectra. 

Recent model calculations by Aboudaram et aL (1984) 
show that the white light flare does not arise from heating 
at T5000 > 1 » because the predicted AI/I^ would then vary 
as 1/X, in contrast to the observations in Figure 3.23. 
Aboudaram et al . , also show that the continuum cannot be 
explained by heating of the upper photosphere and tempera- 
ture minimum region; such models produce AI/I^ increas- 
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ing with wavelength in the visible and, in particular, they 
do not reproduce the observed increase in contrast below 
4000 A. 

Although the spectral data described here are fragmented 
and sometimes ambiguous, we note that the brightest white 
light flares observed with spectrographs (e.g., Donati-Falchi 
et ai, 1984; Neidig and Wiborg 1984), as well as as the 
events observed with filters at several wavelengths, show 
strong Balmer continuum which must arise in a 
chromospheric-like regime. 

In order to estimate the temperature in white light flares 
it is first required to know how the opacity varies with 
wavelength. But even if the latter is known, the interpreta- 
tion will be limited by our incomplete knowledge of the op- 
tical thickness, departures from LIE, and by possible 
systematic errors in the data. If we assume emission, 
where the opacity varies as X^, the observed slope and in- 
tensity of the continuum in the visible region can satisfac- 
torily be fitted to chromospheric temperatures and densities 
by making appropriate (although arbitrary) choices for geo- 
metric thickness and departures from LTE. In this way the 
white light flare continuum observed by Machado and Rust 
(1974) yielded temperatures of 8500-20,000 K. Similarly, 
the April 24 flare (well-observed with broad-band filters) can 
be fitted to 10,(XX)-20,(XX) K. On the other hand, if H" con- 
tributes to the white light flare continuum, then we might 
expect Tg < 8000 K for the layers in which this emission 
arises. 

The temperature estimates discussed here do not allow 
for effects due to possible unknown sources of opacity in 
the blue. If such opacity is present, no further progress can 
be made until these sources are identified. Finally, we again 
remind the reader that the spectrograms obtained thus far 
may not accurately represent the brightest kernels of white 
light flares; hence a precise description of the brightest white 
light flare spectra, including possible contribution by line 
emission, remains in question. 

3.2.6.3 Timing Relationships 

In Figure 3.25 we compare the hard X-ray (HXR) burst 
profiles with the total power radiated in the white-light con- 
tinuum for six white light flares. One event (Feb. 1 1 , 1980) 
was obtained from Big Bear Solar Observatory (BBSO) at 
a single wavelength (4306 A); one event (July 1, 1980) is 
based on data from both BBSO (X 3862 A) and NSO/SP 
(XX 3610, 4275, 4957, 5650 and 6203 A); the remaining 
four events use NSO/SP data at five wavelengths. In com- 
puting the total power it is assumed that the white light flare 
emission extends over the range 2500-10,000 A . Time reso- 
lution varies from 15 s to 30 s except for the July 1 event 
(Zirin and Neidig 1981) where the resolution was 90 s in 
the five wavelengths of the NSO/SP data and 15 s in the 
BBSO data. In the case of the February 1 1 event we assume 


that the white light flare spectrum is qualitatively similar to 
the spectra in Figure 3.23. 

The most striking feature of Figure 3.25 is the overall 
similarity between the white light and HXR. With the ex- 
ception of the late phase of the July 1 flare, the correlations 
suggest that the HXR and white-light sources may be closely 
related, especially in the impulsive phase. This result con- 
firms the conclusions of an earlier study (Rust and Hegwer 
1975) which compared HXR and white-light intensity data 
at two wavelengths (4950 and 5900 A). 

The most unfavorable of the six comparisons in Figure 
3.25 is the July 1 flare. In this case there is good agreement 
between hard X-rays and white light in the impulsive phase, 
but the brightest white light flare point occurred later and 
approximately coincided with the peak thermal energy con- 
tent (Antonucci 1981) and Fe XXV emission (Phillips 1981). 
This suggests an association with the hot (T » 20x 10® K) 
thermal plasma in the flare. Similarly, Zirin and Tanaka 
(1973) observed a white light flare wave event in the hot ther- 
mal phase of the August 2, 1972 flare. 

In summary, the HXR and white light are well-correlated 
in most cases, while in others the white light rises with the 
HXR and declines with the hot thermal plasma. In this sense 
the behavior of the white light in the impulsive and thermal 
phases bears some similarity to the rise and fall of Ha emis- 
sion relative to hard and soft X-rays in non-white light flares 
(Zirin 1978). 

3.2.6.4 Energetics 

The principal datum that motivates our study of white 
light flares is the large observed continuum flux. For nine 
flares the mean peak flux is estimated to be 1.5 X 10*® ergs 
cm"2 s“* between 2500 and 10,000 A (Neidig and Cliver 
1983), with the largest reported flux for a single event be- 
ing 5 X 10*® erg s“* cm"^ (Slonim and Korobova 1975). Al- 
lowing for the low spatial resolution of these observations, 
we suggest, for purposes of model calculations, that a 

I.. 1 1 . — 1 

A uppvi lAiiiii A/11 UA^ lAUA IS lU” ^^111 - S *. 

The peak continuum flux far exceeds the combined flux of 
all the Balmer lines at the same location in the flare, although 
the line emission, being more widespread, may be compar- 
able in total energy when integrated over time and space. 
The total peak continuum power, obtained by integrating over 
the white light flare area, is less sensitive to resolution, and 
is approximately 10^® erg s"* for the largest white light 
flares. 

These numbers are quite impressive, and they place 
severe constraints on the energy transport mechanisms in 
white light flares. A number of transport processes have been 
suggested in the literature, although none can be claimed to 
be satisfactory in general. Chromospheric heating by fast 
electrons has been a long-time favorite mechanism of model- 
ers, and we shall demonstrate here that this mechanism might 
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Figure 3.25 Comparison of total white-light power and hard X-ray profiles for 
six white light flares. Optical data are derived from broad-band filtergrams, as 
described in the text. 


work in the impulsive phases of two well-observed white light 
flares, but fails in the post-impulsive (hot thermal) phase of 
one of these flares. 

For both the impulsive and post-impulsive phases of the 
July 1 and April 24 flares, we list in Table 3.3: (1) the ob- 
served peak power in white light, (2) the electron energy 
Eq above which the power in the nonthermal electron spec- 
trum (in a thick target approximation) is equal to the white- 
light power, (3) the stopping height h(E^ of an electron with 
energy in the non-flaring atmosphere, (4) the column den- 
sity nz(E^) of ionized gas that can be penetrated by an elec- 
tron with energy E^, and (5) the estimated column density 
nz(WL) required to produce the observed white-light inten- 
sity via emission at an assumed temperature of 10^ K 
and a density equal to the density in the non-flaring at- 
mosphere at the electron stopping height. The nonthermal 
electron data are derived from Dennis (1981), Batchelor 


(1984), Kane (1983), and Kane et al (1985). We see that 
in the impulsive phase of both flares sufficient power in non- 
thermal electrons might, within limits of observational un- 
certainties, be transported through a column comparable in 
extent to that required to produce the amount of white-light 
emission. The only problem is producing the observed white 
light flare spectrum. Sufficient power is marginally avail- 
able also in the post-impulsive phase of the 24 April flare, 
although for the post-impulsive phase of the July 1 flare the 
nonthermal electron model seems to fail completely. 

Heating by E ^ 1 Mev protons in the impulsive phase 
white light flare was suggested by Najita and Orrall (1970) 
and Svestka (1970). The protons certainly have enough 
range, and the calculations of Lin and Hudson (1976) show 
that, at least in the August 4, 1972, flare, the proton flux 
alone was probably adequate to power a typical white light 
flare at column depths up to 10^'; cm"^. More recent work 
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Table 3.3 White Light Flare and Energetic Electron Transport Parameters 


July 1, 1980 WL Power 

Eo 

h(Eo) 

nz(Eo) 

nz(WL) 

Apr. 24, 1981 WL Power 

Eo 

h(Eo) 

nz(Eo) 

nz(WL) 


by Hudson and Dwivedi (1982) has shown that sufficient 
heating by protons may be possible even in the upper photo- 
sphere (column depth 5 x 10^2 cm" 2 ), but that the heating at 
depths near rsooo = 1 (column depth 2 xlQ 2 ^ cm" 2 ) is 
negligible. Lin and Hudson originally dismissed the proton 
heating mechanism because of an apparent delay in die 2.2 
MeV 7 -ray lines relative to the flare impulsive phase. But 
recent SMM data shows no delay in the energetic protons 
that produce (via collisions) the free neutrons responsible for 
the 2.2 MeV 7 -ray line, so the delay is not a problem. This 
conclusion is corroborated by Ryan et aL (1983) who show, 
in addition, that while proton heating in the impulsive phase 
of the July 1 flare might have been possible, the proton event 
had ended before the post-impulsive white-light emission 
reached maximum. Thus, we conclude that heating by non- 
thermal electrons or protons cannot account for the post- 
impulsive phase white light flare. 

A more recently investigated mechanism for the white 
light flare is a two-stage transport process. We know that 
sufficient energy is available in E ^ 35 keV electrons even 
in the post impulsive phase. These electrons are stopped near 
the 1000 km level, and the problem then is to efficiently trans- 
port the energy to regions of higher density to produce the 
optical continuum. Livshitz etal (1981) have concluded that 
a shock wave may propagate downward from the heated 
chromospheric region with velocity 1000 km s"*, and that 
this could produce heating at lower levels where the white 
light flare is seen. The principal difficulty here is the rela- 
tively short duration of the shock-induced effect 
(^ 10 s), which is in contrast to observed white light flare 
durations as long as 10 minutes for a single bright kernel 
[see Figure 3.25, and Kane et al (1985)]. 

The post-impulsive white light flare emission in the July 
1 flare seems to be associated with the hot thermal phase 
of the event. We might be led, therefore, to consider heat 
conduction as a transport mechanism, were it not for the fact 


Impulsive Phase 

Post-Impulsive 
(Thermal Phase) 

2.3 X 102’ 

ergs s"^ 

4.5 X 102’ ej.g g- 

53 keV 


35 keV 

880 km 


1020 km 

1.2 X 102* 

cm“2 

5.2 X 102® cm-2 

3.9 X 102* 

cm‘2 

2.1 X 1022 cm-2 

1.4 X 102’ 

ergs s‘^ 

7.0 X 102’ g-i 

119 keV 


85 keV 

650 km 


730 km 

6.0 X 102* 

cm“2 

3.1 X 102* cm-2 

3.3 X 102® cm-2 

3.3 X 102* cm-2 


that sufficient conductive flux can be obtained only at tem- 
peratures of millions of degrees in the white light flare source. 
Clearly, this is incompatible with a chromospheric origin for 
the white light flare; moreover, it requires an emission meas- 
ure several orders of magnitude larger than is observed from 
soft X-ray data. Thus we are unable to identify a satisfac- 
tory mechanism for the energy transport in post-impulsive 
phase white light flares. 

3,2.7 Ha Emission 

Study of the morphology, timing and spectrum of Ha 
(T 1Q4 K) emission during the impulsive phase, relative 
to hard X-rays and microwaves, bears on several basic ques- 
tions of flare energy transport. Do flare heating effects oc- 
cur nearly simultaneously (say within 1 s) throughout loops, 
down to Ha formation depths? Do fast electrons penetrate 
as far as the chromosphere? Is the observed depth depend- 
ence of flare heating in the chromosphere the dependence 
predicted by stopping of nonthermal electrons? Since the 
Skylab era much progress has been made at both observa- 
tional and theoretical approaches to these questions. 

3.2.7. 1 Observations 

During the recent solar maximum the usefulness of Ha 
observations was increased markedly by improvements in 
observing techniques that ensured, for the first time, during 
the impulsive phase, observations with high temporal, spatial 
and spectral resolution, coordinated with temporally, 
spatially, and spectrally resolved X-ray and microwave data. 

a. Timing, The temporal resolution and precision of 
measurement in Ha and either hard X-rays or microwaves 
required to discriminate between conduction front and ener- 
getic particle models of the transport of energy from the co- 
rona to the chromosphere during the impulsive phase is 
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approximately 1 s. It is also necessary to have spatial reso- 
lution, since not all parts of the chromospheric flare neces- 
sarily reflect the same transport process. In the past, the lack 
of adequate temporal resolution led to much confusion about 
the time delays between hard X-rays (or microwaves) and 
Ha (cf. Vorpahl 1972, Zirin 1978). 

The first observations with sufficient resolution to be 
compelling have been obtained recently by Kaempfer and 
Schoechlin (1982) and Kaempfer and Magun (1983). In a 
study of one flare with 1.4 s and 100 ms temporal resolu- 
tion in Ha and microwaves respectively, Kaempfer and 
Magun found evidence for both fast electron transport y at 
one site of a flare, and hydrodynamic or nonclassical con- 
ductive transport^ at other sites of the same flare. At the 
former site, they observed Ha and microwave synchronism 
to within two seconds; the lack of delays at different micro- 
wave frequencies also supported an energetic electron in- 
terpretation. At the latter sites, delays of about 10 s were 
observed. The authors show that these delays are consistent 
with the propagation of disturbances at about 2000 km s‘^ 
i.e. roughly the velocities expected for collisionless conduc- 
tion fronts (Brown, Melrose and Spicer 1979). 

b. Morphology. The spatial coincidence of Ha and hard 
X-ray emission during the impulsive phase of Type B flares 
supports the idea that both emissions often arise as a conse- 
quence of the transport of coronally accelerated nonthermal 
electrons to the footpoints of flare loops (see Section 3.2.2). 
Although this conclusion applies only to impulsive hard X-ray 
flares, not gradual ones, there is no controversy over this 
spatial relationship in the horizontal plane. 

On the other hand, there has been a continuing con- 
troversy over the height of origin of impulsive chromospheric 
flare emission. On one side, Zirin (1978) has argued that 
impulsive Ha emission is produced by direct collisional ioni- 
zation and recombination in an elevated source— prominence 
material— well above the quiet chromosphere. This belief is 
based on line center Ha filtergrams at and near the limb, 
during hard X-ray bursts. Zirin ’s observations show several 
events in which impulsive phase Ha brightening is seen in 
structures that typically extend 10^ km above the limb. On 
the other side, theoretical models all show that the Ha emis- 
sion of flares comes about primarily thermally, from deeper 
than the quiet chromosphere (see, e.g., Brown, Canfield and 
Robertson 1978 or Canfield, Gunkler and Ricchiazzi 1984). 
The latter picture is also what is inferred for the primarily 
chromospheric source of impulsive 10- 1030 A EUV emis- 
sion (Donnelly and Kane 1978), and is discussed in Section 
3.2.5). 

A recent observation of a limb flare, by Kurokawa 
(1983), bears on the controversy. It shows that flare emis- 
sion sites observed at Ha line center ± 2.4 A in a limb flare 
show abrupt intensity changes directly correlated with micro- 
wave fluctuations. It also shows that these emission sites are 


brightest at Ha ± 1 .2 A , where they appear about 18(X) km 
above the photosphere, but are obscured by the quiet chro- 
mosphere at Ha line center. Kurokawa (1983) concluded that 
the main part of the flare chromosphere is confined to the 
low chromosphere, much lower than the surrounding un- 
disturbed structures. 

c. Spectra. Whereas before the recent solar maximum 
there were observations of spectral line profiles from flare 
kernels during the impulsive phase of perhaps one flare (Zirin 
and Tanaka 1973), there are now observations of hundreds. 
Many of these are accompanied by X-ray and microwave 
observations. Such observations of the Ha line (Acton et al , 
1982, Gunkler et al , 1984, Canfield and Gunkler 1984) pro- 
vide strong evidence for the propagation of nonthermal elec- 
trons down into the chromosphere during the impulsive 
phase. 

Figure 3.26, from Canfield and Gunkler (1984), shows 
Ha line profile evidence for chromospheric heating by non- 
thermal electrons during the flare of 7 May 1980, Two of 
the times shown, 14:56: 16 and 14:56:42 UT, are during the 
impulsive phase, and the hard X-ray emission at all ener- 
gies above 30 keV was 3-4 times stronger at the second time. 
There were two distinct kernels in this flare, shown by the 
cross-hatching in Figure 3,26; 2-3 times more 16-30 keV 
emission was observed to come from the south (lower) ker- 
nel, During the time of impulsive hard X-ray emission, strong 
Ha wings develop. They disappear very quickly after the 
impulsive emission ends. The number of pixels that show 
such behavior rises and falls in direct proportion to the in- 
tensity of the hard X-ray emission. 

3.2.T.2 Interpretation of Spectra 

Since fast electrons lose several orders of magnitude more 
energy to heating than to bremsstrahlung X-ray photon 
production, one would expect to be able to see clear effects 
in thermally generated chromospheric emission. Ricchiazzi 
and Canfield (1983) recently developed static models of the 
effects of enhanced fluxes of nonthermal electrons, treated 
strictly collisionally, as well as thermal conduction and en- 
hanced pressure from the flare corona, on the temperature, 
density, and ionization structure of the chromosphere. Their 
work assumed that the electrons last long enough to estab- 
lish a new hydrostatic equilibrium in a closed loop (at least 
tens of seconds). Subsequent work, in an impulsive approx- 
imation, by Canfield, Gunkler and Ricchiazzi (1984), hen- 
ceforth CGR, assumed that very little time (at most a few 
seconds) has passed since the start of electron heating, so 
no mass motions have had time to develop. These models 
superseded work carried out in the Skylab era by Brown, 
Canfield and Robertson (1978). 

By using these model flare chromospheres to compute 
profiles of optically thick spectral lines, which are formed 


14:56:42 UT 


14:57:07 UT 


14:57:33 UT 


Figure 3.26 Impulsive-phase Ha profile observations. The two large contiguous squares framing each of the 
6 panels indicate two 8” x 8' HXIS pixels; the smaller rectangles indicate 2" x 2.67" Ha pixels. The shaded 
Ha pixels indicate the north and south kernels. In each Ha pixel its spectrum at the indicated time l(AA) is plotted 
in units of 1^, the observed quiet sun continuum intensity near the flare site. Each pixel's preimpulsive phase 
spectrum (145525 UT) is shown dotted, for comparison. From Canfield and Gunkler (1984). 
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over a range of depths along the direction of propagation 
of the electrons, one can predict theoretical line profile sig- 
natures of the distribution of heating appropriate to fast elec- 
trons. CGR calculated Ha line profiles based on both 
hydrostatic and impulsive model chromospheres. They 
showed that both hydrostatic and impulsive model at- 
mospheres heated by nonthermal (power-law) electrons 
characteristically produce wide and bright Ha profiles with 
a central reversal. Enhanced thermal conduction reduces the 
width and total intensity of the profiles. High thermal con- 
duction alone cannot account for flare Ha enhancements. 
High coronal pressure dramatically increases the width and 
total intensity of the Ha profiles, while reducing the central 
reversal. 

CGR found that only high values of the flux of nonthermal 
electrons produced Ha profiles with obvious broad (Stark) 
wings. Their results are shown in Figure 3.27. In the upper 
panel, Ha profiles are shown for the impulsive model chro- 
mospheres; those in the lower panel are for hydrostatic 
models. In both panels it can be seen that pronounced wings 
develop only for values of electron energy flux (above 20 
keV) above about 10*® ergs cm'^ s“*. The calculations in- 
dicate that no combination of values of thermal conduction 
and coronal pressure can give the same strong Stark wings, 
accompanied by central reversal, as those associated with 
nonthermal electron heating. 

The observations cited above show that Ha profiles de- 
velop obvious broad wings only in spatial and temporal coin- 
cidence with hard X-ray emission. Gunkler etal (1984) and 
Canfield and Gunkler (1984) found, by quantitatively com- 
paring the extent of observed Ha wings to the theoretical 
profiles of Figure 3.27, that values of F 20 10** ergs 

cm"2 s“* were implied. They showed that this same value 
was implied by combining the power of electrons above 20 
keV inferred from the HXRBS observations with the ob- 
served Ha stark-wing area. From the point of view of trans- 
port theory, the implication of this result is that the electrons 
that produce impulsive hard X-rays indeed have a range that 
is in at least order-of-magnitude agreement with the predic- 
tions of a strictly collisional thick-target model of electron 
propagation. The fact that this inferred value of F 20 is so large 
(approaching the return-current stability limit) is interesting. 
This is not the first time such large electron fluxes have been 
found from chromospheric spectra; using an inversion 
method, F 20 — 10*^ ergs cm"^ s"* was inferred from semi- 
empirical flare chromospheric models (Machado et ai , 1980) 
by Emslie, Brown and Machado (1981). 

3.3 THEORETICAL STUDIES OF TRANSPORT 
PROCESSES 

3.3.1 Electron Beams and Reverse Currents 

The existence of large numbers of nonthermal electrons 
in solar flares has been supposed for some time now as means 


of explaining, through the mechanism of electron-proton 
bremsstrahlung, the high flux of hard X-rays observed in 
large events (see e.g., Brown 1971). With observations 
(Tanaka, Nitta, and Watanabe 1984) now showing the hard 
X-ray continuum (i.e. line-subtracted) spectrum extending 
in an unbroken power-law down to energies ^ 7 keV, it 
appears that, at least in some events, suprathermal 
bremsstrahlung-producing electrons exist at all energies 
above this value. This implies that an extremely large num- 
ber of these electrons must be accelerated in the primary 
release process, whether through directed acceleration (“non- 
thermal” or “beamed” model) or bulk energization (“ther- 
mal” model). In particular, in the nonthermal model such 
hard X-ray observations lead to nonthermal electron injec- 
tion rates > 10^^ s"* above 25 keV (Hoyng, Brown and 
van Beek 1976) and so, by extrapolation of typically observed 
power-laws, > 10^® above - 7 keV. 

One of the major theoretical problems of the last decade 
has been to describe quantitatively the energy, momentum 
and charge transport effected by such strong beams. Early 
treatments (e.g., Brown 1972, 1973; Emslie 1978) neglected 
collective effects within the electrons comprising the beam, 
and instead concentrated on the mean behavior of individual 
electrons through Coulomb interactions with the ambient 
plasma particles. It was soon recognized, however, (Hoyng, 
Brown and van Beek 1976; Knight and Sturrock 1977; 
Hoyng, Knight and Spicer 1978) that the passage of such 
a large number of electrons through the stationary solar at- 
mosphere constituted an extremely large current. This re- 
quired that a reverse (or return) current be set up in the 
ambient plasma, in order to achieve both charge and cur- 
rent neutrality, and to provide fresh electrons to the acceler- 
ation site. We refer the reader to Knight and Sturrock (1977) 
for a discussion of the unacceptably large electrostatic and 
magnetic fields produced in the absence of a reverse cur- 
rent, and to Hoyng, Brown and van Beek (1976) for a dis- 
cussion of the “electron number” problem associated with 
the depletion of electrons from the coronal acceleration site. 
This return current must be driven through a medium of finite 
resistivity, modifying (by an as yet unresolved amount— see 
below) the energetics of the primary electron beam. In ad- 
dition, other collective effects, such as two-stream instabili- 
ties and collisionless pitch-angle dispersion, modify the 
dynamics of such beams, with implications for both their hard 
X-ray production and atmospheric heating. Finally, the beam 
may have yet another significant consequence beyond those 
of photon production and heating: the stopping of the beam 
may effect significant acceleration of the coronal and chro- 
mospheric plasma through which it passes. This section dis- 
cusses the following aspects of nonthermal electron beams: 
the energetics of, and driving mechanism for, the return cur- 
rent, and collective instabilities driven by the beam itself. 
The motivation for such in-depth study of nonthermal elec- 
tron beams is provided by recent observations from SMM 
and other spacecraft, which indicate a substantial nonther- 
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EFFECT OF NONTHERMAL ELECTRONS 



EFFECT OF NONTHERMAL ELECTRONS 
(NOMINAL CORONAL PRESSURE) 



Figure 3.27 The effect on Ha line profiles of varying the energy 
flux of nonthermal electrons above 20 keV, F 20 . The electron beam 
is assumed to be vertical, and the spectral index of the power-law 
electron energy distribution function is 5 = 5. Intensity is meas- 
ured in units of the quiet sun preflare continuum near Ha. Fs is the 
conductive flux from the corona, measured at T = 10* K. Po is the 
gas pressure at the top of the coronal loop, the assumed accelera- 
tion site. Top: impulsive model. Bottom: hydrostatic model. From 
Canfield, Gunkler and Ricchiazzi (1984). 






mal electron injection in some events (Brown, Hayward and 
Spicer 1981; Hoyng 1981; Hoyng et al, 1981). 

3.3. 1.1 Energetics of the Reverse Current 

Emslie (1980) showed that the ohmic energy losses sus- 
tained in driving the reverse current through the finite 
resistivity of the ambient plasma could, in some cases, ex- 
ceed the losses suffered through direct collisions of the beam 
electrons with ambient ones. For classical resistivity, the ratio 
of ohmic to collisional heating scales like F n"^ 
where F is the injected electron flux, n the ambient density 
and Tg the electron temperature (Emslie 1980), and hence 
only for relatively large fluxes of electrons passing through 
relatively cool coronae could this ohmic heating dominate. 
This raises the question of the self-consistency of Emslie’s 
(1980) treatment through considerations of energy balance 
and, perhaps more fundamentally, through the requirement 
that the electron drift velocity associated with the reverse 
current be below the threshold for plasma instabilities to de- 
velop, typically around the ion-sound speed c^ = 
(kgTg/mp)*^ 2 ^ where kg is Boltzmann’s constant and the 
proton mass (Fried and Gould 1961). These matters were 
addressed independently by Emslie (1981) and by Brown and 
Hayward (1982). Emslie (1981) considered the flux limita- 
tion imposed by the stability requirement on the reverse cur- 
rent drift velocity. He showed that this could reduce the 
importance of reverse current ohmic heating, except during 
the very early stages of the event, when low preflare coronal 
temperatures result in a very large plasma resistivity (see also 
Knight and Sturrock 1977), and also during the late phase 
of long events (t > 10 s), by which time the ambient pro- 
tons could be significantly heated by the hot electrons, result- 
ing in an equilibration of ion and electron temperatures, an 
increased threshold drift velocity (Fried and Gould 1961) and 
so a larger possible ohmic to collisional energy loss ratio. 
Brown and Hayward (1982, see also Hayward 1984) con- 
sidered the other aspect of the self-consistency problem, 
namely the determination of the ambient plasma electron tem- 
perature Tg through consideration of the energy equation. 
Including only dominant terms, this can be expressed as 

r;j2 + V • Fg = 0, (3.6) 

viz. an energy balance between ohmic heating and thermal 
conductive cooling. Here r/ is the plasma resistivity, j the 
beam (or reverse) current density, and Fg the thermal con- 
ductive flux. Equation (3.6) assumes that direct collisional 
losses are unimportant, an assumption that can be tested a 
posteriori by comparing the ohmic losses in a plasma whose 
temperature follows from Equation (3.6) with the collisional 
losses in the same plasma. The collisional losses in turn de- 
pend on the density n and the average individual electron 
energy E. Thus the ratio of ohmic energy losses to collisional 
ones depends on F, E, Tg and n. In addition, from the values 


of n and Tg [found from solution of Equation (3.6)], we can 
test whether the beam flux exceeds the stability threshold 
nVgrit’ where depends solely on Tg. From these consider- 
ations, Brown and Hayward (1982) deduced that when the 
reverse current is stable to the generation of plasma tur- 
bulence it is unimportant energetically; the effects of an un- 
stable reverse current are yet to be determined. 

Emslie (1985) has recently contested the Brown and Hay- 
ward (1982) results, on the basis of the inadequacy of their 
treatment of the conduction term in (3.6) as 
where is the Spitzer (1962) coefficient of thermal con- 
ductivity and Lg 3 the stopping distance of the beam electrons 
under the decelerating force produced by the charge separa- 
tion electric field E = r;j. He argued that the boundary con- 
dition determining the point at which the electrons stop is 
not in fact Lg^, but instead the half-length L of the flare loop 
(beyond which Coulomb collisions in the high density chro- 
mosphere effectively stop the beam in a very short distance), 
and therefore recast Equation (3.6) in the form 

rjoT' -h d/dz {xjp dTg/dz) = 0, (3.7) 

imposing a boundary condition on Tg at a fixed distance L 
corresp>onding to the half-length of the flare loop. Using this 
revised version of the energy equation, and also a more realis- 
tic power-law-type injected electron energy distribution, 
Emslie (1985) found, in contrast to the results of Brown and 
Hayward (1982), that a regime (albeit a small one in the avail- 
able parameter space) where a stable reverse current 
dominates the energetics of the upper coronal plasma can 
exist. Physically, this is because the shorter and more phys- 
ical boundary values of L imposed by Emslie result in a much 
lower apex temperature; since the resistive losses for a stable 
reverse current scale as Tg~3/2 ^nd the drift velocity stabil- 
ity threshold scales as Tg*^^^ a reduction in this peak tem- 
perature can have relatively large effect on the importance 
of ohmic heating, while having a relatively small effect on 
the maximum beam flux allowed to pass stably. For low den- 
sities, this stable reverse current heating is dominant, while 
for high densities collisional heating dominates. Thus there 
is a minimum energy input e = max (ecoiiisions* ^reverse cur- 
rent) into the coronal plasma, and, in turn, a minimum peak 
coronal temperature for a given hard X-ray flux, independent 
of coronal density. This relationship between hard X-ray 
burst intensity and peak coronal temperature (deduced from 
soft X-ray observations) in different events may afford ob- 
servational tests of (a) the thick-target model and (b) the rela- 
tive role of reverse current and collisional heating in such 
a model. 

3.3. 1.2 Driving Mechanism 

A fundamental issue for the physics of nonthermal elec- 
tron beams is whether the reverse current is established elec- 
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trostatically or inductively. Spicer (1982) and Spicer and 
Sudan (1984) have claimed that an assumption made by all 
previous authors, namely that the reverse current is in steady- 
state balance with the injected beam flux, and that it is es- 
tablished by the electrostatic field associated with the charge 
separation effected by the beam, is incorrect. Rather, they 
argued that the unneutralized current of the beam sets up, 
through Ampere’s law, a self-magnetic field B which then 
establishes a return current inductively through Lenz’s law. 
In addition, they argued that this return current decays resis- 
tively in the ambient plasma, so that eventually it can no 
longer effectively neutralize the beam; the resulting “bare 
beam” cannot propagate (for the reasons discussed by Knight 
and Sturrock 1977) and so turns off. They suggested that 
this may be the cause of the extremely rapid decays found 
in some hard X-rays events (Kiplinger et al , 1983). Brown 
and Bingham (1984) have criticized Spicer and Sudan’s work 
on the grounds that they neglected the displacement current 
term in Ampere’s law at the head of the beam. The inclusion 
of this term, according to Brown and Bingham, returns us 
to a scenario where the reverse current is indeed set up elec- 
trostatically, over a Debye length scale at the head of the 
beam; since the time required for the beam to traverse such 
a small distance is much smaller than timescales for hard 
X-ray bursts, they argued that a steady-state description 
(Knight and Sturrock 1977; Emslie 1980) is indeed valid. 
With regard to the resistive decay of the reverse current, they 
agreed with Spicer and Sudan that the finite curl of the E 
field (due to edge effects at the outside of the beam) leads 
to the establishment of an azimuthal magnetic field B_l. While 
this field would indeed choke off the beam over a short time- 
scale (Alfven 1939; Lawson 1957) in the absence of a strong 
guide field, Brown and Bingham (1984) pointed out that the 
presence of such a guide field, as exists in flare loops, results 
only in a increased helicity in the guiding field lines; as long 
as the azimuthal B field remains small or comparable to the 
guiding longitudinal field Bj|, then little modification to the 
beam energetics results. The consensus reached is that, on 
the one hand, for beams spread out over large areas, edge 
effects are of minor importance, and so the predominant 
mechanism controlling the reverse current is the electrostatic 
field established by the displacement current at the head of 
the beam. On the other hand, for beams of small area, e.g., 
if the injection occurred over many sub-resolution elements 
(a “pepper-pot” scenario), the inductive edge effects would 
predominate. In such a case, “choking-off ’ of the beam 
(Spicer and Sudan 1984) may occur on timescales relevant 
to observations, but for this effect to have an appreciable 
hard X-ray signature, many elements of the “pepper-pot” 
would have to act coherently, and it is far from clear why 
this should occur. However, in both these cases it appears 
that the earlier steady-state electrostatic treatments of the 
reverse current energetics give valid results. 


3.3. 1.3 Collective Instabilities 

Emslie and Smith (1984) showed how, due to the inverse 
square energy dependence of the Coulomb collision cross- 
section, the velocity distribution for a nonthermal electron 
beam injected into a cold plasma is preferentially depleted 
at low energies, resulting in the formation of a “hump” in 
the combined (background plus beam) distribution function. 
This type of distribution is well-known to be unstable to the 
growth of Langmuir waves; Emslie and Smith calculated the 
resulting wave level and discussed nonlinear interactions of 
these waves, through the process f + — t, (where f is a 
Langmuir wave and t a transverse wave with angular fre- 
quency » 2o)pe, where Wpe is the electron plasma frequen- 
cy), to produce microwave radiation. They showed that the 
level of microwave radiation produced by this process is, 
for typical thick target beam parameters, much larger than 
observational upper limits on microwave fluxes from flares. 
This implies that either the nonthermal microwave radiation 
produced is strongly absorbed by gyroresonance absorption 
in the surrounding plasma, with corresponding strong im- 
plications for the energetics of the flare as a whole (Melrose 
and Dulk 1982b), or that the thick target electron heated 
model in its present form is incorrect. 

Mok (1985) has examined another instability connected 
with the scattering of electron beams in the solar atmosphere 
— the anomalous Doppler resonance instability (Kadomstev 
and Pogutse 1968). This instability is driven by anisotropy 
in the electron velocity distribution f(V||,NQ^) and operates 
through “beat” resonances of the form k*v = co + nQ, 
where k is the wavenumber, v the electron velocity, w the 
electron plasma frequency, the gyrofrequency and n an 
integer. Using a power-law pitch angle distribution centered 
in the direction of the guiding magnetic field lines, with a 
moderate degree of anisotropy (0 * 0.25), Mok found that 
for low Wpe/n (e.g., in the corona), no resonant amplifica- 
tion of waves through this instability results. However, with 
the increase in cope caused by increasing density in the chro- 
mospheric regions of the loop, the ratio o)^/Q becomes large 
enough so that the instability threshold is exceeded. This 
results in a growth of lower hybrid waves (Krall and Trivel- 
piece 1973), leading to fast collisionless scattering of the 
beam electrons into a more isotropic distribution. This has 
several effects. First, it causes the beam to suddenly slow 
down. Second, there is strong local heating produced near 
the region where the onset of instability occurs, possibly af- 
fecting the ratio of hard X-ray and EU V fluxes (cf . Section 
3.2.5). Third, some microwaves are emitted in this region, 
although they may not represent a dominant contribution to 
the microwave emission during the impulsive phase. 

3.3. 1.4 Electron Beam Momentum 

The strong observational evidence for the applicability 
of the thick target electron beam model to the transport of 
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impulsively accelerated electrons has led many authors to 
carry out studies of its hydrodynamic and energetic conse- 
quences on the structure and dynamics of the target at- 
mosphere (see, e.g., Sections 3.2.5 and 3.2.7). However, 
these studies have all focused their attention on the heating 
aspects of the electron beam; they have neglected the effect 
of direct momentum deposition in the equation of motion. 
As pointed out by Brown and Craig (1984), beam momen- 
tum is potentially of great interest to impulsive phase trans- 
port because it is transferred directly to the atmosphere, and 
hence produces an immediate acceleration ag, whereas ac- 
celeration by nonhydrostatic thermal pressure gradients, ap. 
takes time to build up. 

Brown and Craig found that an electron beam, for high 
values of the energy flux, produces a beam acceleration that 
is much greater that solar gravity (ag*? 100 g^ , for beam 
energy flux 10^* ergs cm"^ s"^. They compared ag to ap 
on the assumption that beam heating is balanced by conduc- 
tion. Their comparison implied that ap was less that ag for 
classical conduction and of the same order as ag for satu- 
rated heat flux. Hence, they concluded that beam momen- 
tum deposition is important even after nonhydrostatic thermal 
pressure gradients have developed. 

McClymont and Canfield (1984) came to quite a differ- 
ent conclusion. They estimated the time scale for ap to ex- 
ceed ag, arguing that conduction and radiation can be 
neglected on the short time scales of interest, and so con- 
sidered beam heating alone. They found that in the corona, 
ag dominates ap for at most the first few seconds; in the 
chromosphere it will do so even more briefly (^ 10“^ s). 
The basic reason for this is the large density gradient across 
the transition region in the preflare atmosphere. Since colli- 
sional heating gives rise to a temperature increase rate which 
is (locally) independent of density (Emslie 1978), the rapid 
heating of the preflare transition region layers to coronal tem- 
peratures gives rise to a large pressure gradient across this 
region. Hence they concluded that electron beam momen- 
tum deposition is not important to the global evolution of 
the flare atmosphere. 

Both Brown and Craig (1984) and McClymont and Can- 
field (1984) also considered the hydrostatic case, i.e. an elec- 
tron beam that has been on for long enough to reach a new 
pressure equilibrium. Because ag can so greatly exceed 
gravity. Brown and Craig concluded that steady-state 
hydrostatic models require revision. However, McClymont 
and Canfield called attention to the fact that flare electrons 
are accelerated in closed loops. Due to containment of the 
heated plasma, the thermal coronal pressure is much greater 
than that due to gravity alone. As a result, the addition of 
beam pressure to gravitational and containment pressure has 
a negligible effect. 

At this time no study has yet been carried out that ex- 
amines the hydrodynamic importance of beam momentum 
deposition, taking fully consistent account of beam heating, 
conduction, radiation, and mass motion, in order to unam- 


biguously demonstrate which of these conflicting views is 
correct. For hydrostatic models, on the other hand, it seems 
clear that the electron beam pressure can safely be neglected. 
However, under rather specialized conditions the momen- 
tum imparted to the atmosphere by proton beams may be 
significant (Brown and Craig 1984, Tamres, Canfield and 
McClymont 1985). 

3.3.2 Proton Transport 

The timing, spectrum and morphology of flare protons, 
when compared to electrons, provide important clues on the 
nature of flare acceleration and transport. Here we concen- 
trate our attention on specific recent advances in understand- 
ing proton transport processes between the acceleration site 
and the lower solar atmospere; the reader interested in in- 
terplanetary proton phenomena is referred to Forman, 
Ramaty and Zweibel (1985). 

3.3.2. 1 Lyman-a Charge-Exchange Emission 

The broad-band X-ray and microwave spectrum of the 
impulsive phase of solar flares has been observed well enough 
so that we can confidently state that the most commonly ac- 
celerated nonthermal electrons are those in the deka-keV 
range. Are the most commonly accelerated protons also in 
this energy range? Or does the fact that gamma rays are the 
only observed emission due to protons mean that flares ac- 
celerate protons primarily to MeV energies? Suppose that 
flares typically accelerate protons, in numbers comparable 
to electrons, in the deka-keV range. What are the observ- 
able manifestations? 

It is known that if energetic protons are injected into the 
largely un-ionized chromosphere, they will efficiently pick 
up electrons from ambient neutral hydrogen atoms (Orrall 
and Zirker 1976). The nonthermal protons thereby become 
nonthermal hydrogen atoms, and radiate the usual hydrogen 
emission spectrum, redshifted by the Doppler effect. Orrall 
and Zirker carried out an equilibrium calculation of charge 
exchange in an empirical model chromosphere that showed 
that the nonthermal La radiation would be readily observ- 
able above the quiet sun Lyman-a spectrum, for an input 
energy flux of 2 x 10^ ergs cm"^ s"^ assuming a power- 
law proton energy distribution with a 10 keV low-energy 
cutoff and exponent 6 = 2.5, 3, and 4. Recognizing that the 
ionization structure of the atmosphere would not long remain 
undisturbed by the energetic protons, they showed that for 
the beam parameters considered, the chromospheric ioniza- 
tion time was in the range 10^ - 10^ s. 

Recently Canfield and Chang (1984) examined the sen- 
sitivity of the Lyman-a photon emission spectrum to proton 
energy. Since the largest flares are thought to inject about 
10“ ergs cm’2 s"^ into the lower atmosphere, they inves- 
tigated the Lyman-a radiation that would be generated by 
this same flux, in the form of monoenergetic proton beams. 
Since the temperature and density structure of the atmosphere 
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is of secondary relevance compared to its ionization struc- 
ture, they assumed a uniform atmosphere. They adopted a 
10% hydrogen ionized fraction, and assumed an equilibrium 
between the fast protons and the ambient chromosphere. 
Their result, shown in Figure 3.28, demonstrates a high level 
of emission over a wide range of proton energies. 

It remains to be demonstrated that the radiation gener- 
ated by charge exchange is produced when the atmosphere 
is heated (and perhaps largely ionized) by the proton beam. 
In particular, for what range of proton beam parameters will 
enough of the atmosphere remain sufficiently neutral to make 
charge exchange emission comparable to that calculated? 
Canfield and Chang (1984) showed that at high input beam 
fluxes the ionization time is much shorter than the charge 
exchange time. At high values of incident proton energy the 
ionization time tends to decrease more rapidly than the charge 
exchange time. At low values a proton accelerated at the top 
of a typical coronal loop does not reach the chromosphere. 
At low values of the incident energy flux F(0) the emission 
predicted by the steady-state calculation is too weak to be 
detectable above a typical preflare (active region) back- 
ground. Finally, at high values the ionization time becomes 


short, making the emission difficult to detect. Notwithstand- 
ing, Canfield and Chang showed that some cases of poten- 
tially great interest (e.g., a beam of 300 keV protons and 
energy flux 10® ergs cm"^ are well approximated by 
an equilibrium calculation for tens of seconds, and easily de- 
tected above the active-region Lyman-a background. Future 
calculations should examine time-dependence, a realistic tar- 
get atmosphere, and spectra including Bessel functions (For- 
man, Ramaty, and Zweibel 1984). 

3.3.2.2 Heating of a Thick-Target Atmosphere 

With the simultaneous operation of the Hard X-ray Burst 
Spectrometer (HXRBS) and Gamma Ray Experiment (GRE) 
on SMM, a few flares exist for which there are good data 
in both hard X-rays and y-rays. These data provide infor- 
mation of the populations of energetic deka-keV electrons 
and Mev protons, respectively. Emslie (1983a) has used the 
available data from both experiments for the 1980 June 7 
flare (see Kiplinger et ai , 1983) to infer the energies in both 
electrons and protons for this event. Through thick target 
modeling of these emissions (Brown 1971, Ramaty 1984), 


NONTHERMAL La SPECTRA FOR MONOENERGETIC PROTONS 



Figure 3.28 Solid curves: theoretical nonthermal Lyman-a spectra for monoenergetic ver- 
tical proton beams of energy flux 10’^ ergs cm"^ s"\ The short-dashed curve, marked 
F, Is the average of two flares at Lyman-a flare maximum (Canfield and VanHoosier, 1 980). 
The long-dashed curve Is an active-region spectrum (Cohen 1982). 
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Emslie has calculated the energy deposited per unit height 
by both electron and proton bombardment, as a function of 
depth in the atmosphere. Apart from the (relatively) unknown 
area factor, determination of the instantaneous energy input 
rate is possible only for electron heating, since hard X-rays 
are “prompt” data, providing information on the instantane- 
ous target-averaged electron flux. On the other hand, some 
7 -rays (e.g., the deuterium formation 2.223 Mev line) are 
produced over a considerable period (compared to observa- 
tional time scales) after the injection of the protons which 
lead ultimately to the 7 -ray emission. We can calculate only 
the event-integrated energy deposition with reliability. 

Proton spectra are in general harder (shallower) than elec- 
tron spectra. Since Mev protons have a greater individual 
energy and stopping depth (Emslie 1978) than deka-keV elec- 
trons, it follows that in general proton heating dominates elec- 
tron heating at large depths, while high in the atmosphere 
electron heating dominates. The “cross-over” point depends 
on the actual shape of the electron and proton spectra, and 
on the relative energy fluxes in the two sets of particles. 
Figure 3.29 shows, for the 1980 June 7 flare, the energy 
deposited per unit height interval by both electron and pro- 
ton bombardment. It can be seen that for this event, elec- 
tron heating dominates all the way down to a particle column 
density N * 3 x cm"^. Since this level is below any 


reasonable location of the chromospheric flare (e.g., 
Machado et aL , 1980; Ricchiazzi and Canfield 1983), it fol- 
lows that electron heating dominates throughout most of the 
region that gives rise to observed flare emissions (with the 
exception of emissions from around the temperature mini- 
mum [see Machado, Emslie, and Brown 1978]). Since the 
1980 June 7 event had a particularly high y-ray to hard X-ray 
flux ratio (Kiplinger 1983), it follows that the latter conclu- 
sion is valid for all flare events (see also Lin and Hudson 
1976), However, Emslie (1983a) pointed out that this con- 
clusion is critically dependent on the assumption of thick tar- 
get emission for both the hard X-rays and the 7 -rays. While 
Ramaty (1984) has conclusively demonstrated that 7 -ray 
emission in solar flares is a thick target process, there is still 
considerable speculation (see, e.g., Brown and Smith 1980; 
Emslie 1983b) on whether hard X-rays are produced prin- 
cipally by thick target bremsstrahlung or by thermal brems- 
strahlung from a confined ensemble of very hot (^ 10® K) 
electrons. In the latter case, the number of precipitating elec- 
trons is reduced relative to the hard X-ray photon flux, and 
less chromospheric heating results (see Emslie and Vlahos 
1980). Thus, if hard X-rays indeed turn out to be emitted 
principally by a hot thermal source, then the importance of 
proton heating in the chromospheric energy budget would 
merit closer examination. 



Figure 3.29 The quantity l^/n = J j J [Q(n,y,z,t)/n(z)l dx dy dt (units ergs per cm of vertical 
height per ambient electron), representing the temporal and areal integral of the specific energy 
deposition rate Q as a function of column depth N = jn(z)dz, for the 1980 June 7 event. 
From Emslie (1983a). 


3-40 



3.3.3 Radiative Energy Transport by 
Amplified Decimetric Waves 

Following the initial work on the application of maser 
theory to radio bursts from the Sun and stars (Holman et al , 
1980; Melrose and Dulk 1982a) and a discussion of the ener- 
getic implications (Melrose and Dulk 1982b), Melrose and 
Dulk (1984) have carried out new work on the transport of 
energy by amplified radio radiation. These transport aspects 
are summarized below; the role of masers in particle acceler- 
ation is discussed in the previous chapter. 

Strong theoretical arguments suggest that up to 50% of 
the energy in a flare should go into radiation at the cyclo- 
tron frequency fig, the exact percentage depending on the 
degree of field line convergence toward the footpoints of the 
flux tubes, and the amount of fast electron energy lost through 
collisions. These arguments (based on an extension of then- 
existing theory by Wu and Lee in 1979) are reinforced by 
consideration of how the electrons accelerated in a flare lose 
the component of their energy perpendicular to the magnetic 
field in order to precipitate into the footpoints of magnetic 
flux tubes and hence heat the chromosphere and generate hard 
X-ray bursts. The perpendicular energy is radiated away 
through amplification of decimetric radio waves. The radi- 
ation travels from the flux tubes of energy release, across 
field lines, and is then reabsorbed at the second harmonic 
of Qg, i.e. at points where the field strength B is half its value 
in the source. From energetic considerations there seems to 
be enough energy in the radio-frequency (RF) radiation to 
heat the plasma of the absorbing region so that it produces 
the soft X-rays. In addition it is likely that this RF heating 
process can account for other properties of soft X-rays in 
the impulsive phases of flares: large source sizes, excess line 
widths, characteristic temperatures, and their temporal 
development. 

3.3.3. 1 Mechanism for Emission and Amplification 

The generation of the high levels of RF radiation is 
thought to occur in an impulsive flare as follows. 

(1) It is assumed that energy release occurs mainly near 
the tops of one or more magnetic flux tubes of the kind il- 
lustrated in Figure 3.30. The accelerated (or heated) elec- 
trons are assumed to have energies E - 10 to 100 keV (or 
T - 10® to 10^ K) and to have a nearly isotropic pitch angle 
distribution, i.e., (^) — (E|j). 

(2) The fast electrons, with v >0.1 c, travel down the 
legs of the flux tubes, those with moderate to large pitch 
angles a. reflecting in the converging fields, and those with 
small a precipitating. In a typical travel time of — 1 s the 
electron pitch angle distribution f(a) in the legs is anisotropic, 
with no upgoing electrons with small a. The anisotropy 
represents a source of free energy because it creates a non- 
equilibrium distribution with an inverted population, i.e., 
there is a region of velocity space where dUdv > 0. 



Figure 3.30 Schematic drawing of a sequence of 
loops showing an energy release region, the maser 
source region in one leg where the field strength is 
B, the cone of maser radiation with reflection of radi- 
ation directed toward higher field strengths, and the 
absorption regions where the field strength is B/2. 
From Melrose and Dulk (1984). 


(3) Provided that the plasma frequency cOp is < Qg, there 
can be a resonant transfer of energy from electrons to waves: 
in particular, electromagnetic waves at the cyclotron fre- 
quency can be amplified. The growth rate of this “cyclo- 
tron maser” is very large, -10^ s~\ so in a few 
microseconds the waves can be amplified by a factor of — e^® 
or so, at which time they saturate, having extracted most of 
the free energy. With B in the expected range of 100 to 1000 
gauss, the amplified radiation is in the range 0.3 to 3 GHz, 
i.e., at decimetric wavelengths. 

(4) Saturation of the maser is expected to occur by 
quasilinear relaxation: the maser converts the perpendicu- 
lar energy of certain electrons into radiation and leaves them 
with mainly parallel energy, i.e. it diffuses them into the loss 
cone at the maximum possible rate. Provided that the den 
sity is low and/or the electron energy is large, so that colli- 
sions are infrequent, this is the only known process which 
can prevent most of the electrons from remaining trapped 
in converging, coronal flux tubes for many bounce periods, 
i.e. for several to many seconds. For example, 30 keV elec- 
trons trapped between mirror points 3 x 10® cm apart would 
bounce - 10^ times (- 10^ s) in a plasma of n^ = 10’ 
cm"^ before collisions would diffuse them into the loss 
cone (Section 3.2.2). 

3.3.3.2 Energy Content 

From studies of hard X-ray bursts it is known that the 
energy going into fast electrons in the impulsive phase is E 
— lO^”^ to 10^® ergs s~^ (e.g., Brown and Smith 1980). 
Provided that the electrons are not highly collimated, about 
half of this is in the parallel component and, again assuming 
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a low collision frequency in the coronal portion of the flux 
tubes, it is deposited in the low atmosphere at the footpoints 
of the loops of energy release. Most of the other half goes 
into the maser radiation, and then (as described below) into 
reabsorption and heating of the coronal plasma in a large 
volume, V — 10^^ to 10^^ cm^, surrounding the loops 
of energy release. Thus the average heating rate is E/V - 1 
to — 10 ergs cm s~\ assuming that the larger values of 
E and V go together. On average this can increase the inter- 
nal energy nkT of the preflare corona by a factor of 2 to 
20 each second. 

3.3.3.3 Reabsorption at Distant Locations 

The RF radiation travels outward from the maser gener- 
ation regions in the manner sketched in Figure 3.30. Only 
one maser emission region is shown, but there are expected 
to be > 10*^ maser pulses per sec, each lasting about 1 ms, 
coming from various points in the legs of the flux tube, and 
from all flux tubes of energy release. Each pulse is narrow- 
band and at a frequency w slightly higher than but in total 
the pulses should cover a frequency range of an octave or 
so because B, hence and co, varies with position in the 
legs of the loops. 

The maser pulses travel across field lines at nearly the 
speed of light until they encounter a region where the field 
strength is half its value at the source: this may be 10^ to 
10"^ km distant, depending on the gradient of B. Then gyro- 
resonance absorption by the ambient electrons occurs at the 
second harmonic, very efficiently, with an optical depth r 
— 10^. A series of pulses can heat a small volume of plasma, 
< 1Q20 cm^, to a few x lO'^ K in - 0.1 s. Temperatures 
above about 3 x 10”^ K energy are difficult to achieve by 
maser heating because of the temperature dependence of the 
collision frequency; since the mean free path is a Tg* ^ 
further energy deposition leads to heating of a larger volume 
of plasma rather producing a higher temperature. 

Turbulent plasma motions should occur because the maser 
pulses heat the ambient electrons in many, small, localized 
volumes, a result of the fine beaming and small absorption 
distance. The hot electrons tend to propagate out of the pock- 
ets of high pressure so created and thus provide an explana- 
tion for the highly broadened soft x-ray lines observed during 
impulsive flares. 

The overall volume heated depends mainly on the field 
gradient in the vicinity of the loops of energy release. The 
model suggests that in a small flare a volume - 10^^ cm^, 
some 10 times larger than the volume of energy release, 
would be heated to ^ 10^ K in 10 seconds. In a large flare 
it is expected that the heated volume would be much greater 
than 10^^ cm^, as would the volume in which energy release 
occurs. Because the maser radiation is initially directed at 
a large angle («70°) to the field lines, most of the energy 
transport is across the field. 


3.3.3.4 Secondary Effects 

According to the RF heating model, the time profile of 
maser emission in simple impulsive flares should mimic the 
hard X-ray profile, and this should follow closely (within 
- 1 s) the profile of energy release. Much of the RF energy 
is deposited in the soft X-ray plasma and is radiated or con- 
ducted away, relatively slowly. Hence the slope of the soft 
X-ray intensity profile is expected to be approximately 
proportional to the hard X-ray intensity as observed. While 
this property is not unique to the RF heating model, it ap- 
plies here to a soft X-ray emitting volume much larger than 
that of initial deka-keV energy release. 

Densification of heated loops and the origin of the blue- 
shifted component of soft X-ray lines can occur in two ways. 
First, energy in the parallel component can “boil off’ chro- 
mospheric material into the loops of energy release. Second, 
as other loops become heated by the maser radiation, en- 
hanced conduction to the footpoints would boil off material 
there, though presumably not as fast. 

The RF heating model has implications for enhancement 
of chromospheric evaporation, a major topic of the follow- 
ing chapter. As pointed out there, chromospheric evapora- 
tion can take place by both direct Coulomb heating from the 
parallel component of the nonthermal electron energy and 
by thermal conduction from energy deposited in the corona 
(possibly a slower process). Hence, if the RF heating 
mechanism heats loops other than those in which direct 
energy release has taken place, chromospheric evaporation 
may take place over a larger area of the chromosphere. 

If the density in a flaring loop increases to the point where 
cop ^ Qg e.g. ng > 10'® cm"^ if B = 300 gauss, the maser 
would stop. Or if the initial density in the flaring loop were 
such that ojp ^ Qg there would be no maser, and probably 
only a gradual phase. However, in this high-density case. 
Coulomb collisions would scatter and degrade the energy of 
all but the highest energy electrons. Thus the conditions 
would be appropriate for a “high coronal flare” without an 
impulsive hard X-ray burst from footpoints. 

3.3.3.5 Conclusions 

As the RF heating radiation is entirely reabsorbed in the 
corona, it is not possible to observe it directly. However, 
there is strong evidence that this maser mechanism is the ex- 
planation for the very intense auroral kilometric radiation 
of Earth and for similar radiation from Jupiter and Saturn. 
From the Sun and stars, direct maser action (or a related 
process involving conversion of amplified waves) produc- 
ing radiation just above the second harmonic is the favored 
explanation for “microwave spike bursts” of very high 
brightness and circular polarization; Enome (1983) has re- 
cently examined some of the properties of such bursts, and 
earlier studies have been made by Droge (1977), Slottje 
(1978) and Zhao (1982). 


3-42 


3.4 SUMMARY 

In this chapter our main interest has been the transport 
of energy, momentum and charge during the dramatically 
nonthermal and energetic impulsive phase. Not surprisingly, 
the technological and operational advantages embodied in the 
SMM and Hinotori instruments, and their collective and 
coordinated use with forefront ground-based optical and radio 
instruments, have brought about many advances in our un- 
derstanding of impulsive phase flare physics. 

The topic most thoroughly explored in this chapter is that 
of the transport of nonthermal electrons. The thick-target 
electron beam model, in which electrons are presumed to 
be accelerated in the corona and typically thermalized primar- 
ily in the chromosphere and photosphere, is supported by 
observations throughout the electromagnetic spectrum. At 
the highest energies, the anisotropy of 7 -ray emission above 
10 MeV clearly indicates that these photons are emitted by 
anisotropically-directed particles. The timing of this high- 
energy 7 -radiation with respect to lower-energy hard 
X-radiation implies that the energetic particles have short life- 
times. For collisional energy loss, this means that they are 
stopped in the chromosphere or below. Stereoscopic (two- 
spacecraft) observations at hard X-ray energies (up to 350 
keV) imply that these lower-energy (but certainly nonther- 
mal) electrons are also stopped deep in the chromosphere. 
Hard X-ray images show that, in spatially resolved flares 
whose radiation consists of impulsive bursts, the impulsive 
phase starts with X-radiation that comes mostly from the foot- 
points of coronal loops whose coronal component is outlined 
by microwaves. Recent X-ray polarization measurements at 
the low end of the deka-keV range show small values, but 
these are still large enough to be compatible with highly col- 
limated particle injection somewhere in the coronal part of 
a loop. Preliminary analyses of combined hard and soft X-ray 
spectra appear to demand the presence of nonthermal elec- 
trons early in the impulsive phase. The thick-target electron- 
beam model accounts for both the close temporal coincidence 
of UV and hard X-ray bursts and the relative decrease of 
UV emission with increasing hard X-ray intensity. White- 
light emission from the largest flares is closely related tem- 
porally and energetically to the nonthermal electrons that 
produce deka-keV X-rays. Moreover, white-light spectra im- 
ply that commonly such emission originates in the chro- 
mosphere, where deka-keV electrons are stopped. Both 
timing and spectra of Ha emission imply that energetic non- 
thermal electrons penetratively heat the chromosphere dur- 
ing the impulsive phase. Finally, we now suspect that 
microwave spike bursts may require a level of microwave 
maser emission that dramatically enhances the precipitation 
of thick-target electrons. 

Although the thick-target nonthermal-electron model 
meets with some success, it describes only a fraction of ob- 
served flares, and shows some puzzling and challenging 
failures. The limited amount of stereoscopic hard X-ray data 


at 350 keV shows little directivity, compared to the predic- 
tions of models of collisional electron transport in loops. If 
this trend is borne out by more comprehensive observations, 
it could provide a strong limitation on the thick-target model. 
The height dependence of the hard X-ray spectrum is sys- 
tematically different from that predicted. The majority of im- 
ages at the low end of the deka-keV energy range show single 
X-ray sources; in many flares, the X-ray producing electrons 
in this energy range are obviously confined to the corona. 
It now seems clear that direct collisions by nonthermal elec- 
trons are not necessary to understand the production of inner- 
shell Fe Ka emission; it appears that impulsive variability 
of these lines comes about because of irradiation of the chro- 
mosphere and photosphere by an impulsively varying spec- 
trum that extends down to at least 7 keV. The thick-target 
model appears to predict much more EUV emission, in 
proportion to hard X-ray emission, than is observed. In some 
flares, white-light emission is observed well beyond the end 
of the impulsive phase, which appears utterly inexplicable 
in the thick-target model; such emission appears to corre- 
late more strongly with the thermal X-ray emission of the 
hot (T > 2x10'^ K) thermal plasma. Furthermore, some 
flares show spectra that imply that the white-light emission 
comes from deeper than the chromosphere, in regions that 
are collisionally inaccessible to deka-keV electrons. Finally, 
it is not well understood whether electron beams of the re- 
quired flux can be maintained, in the presence of return- 
current instabilities known to exist in laboratory plasmas, 
and in view of the observational uncertainty in beam area. 

Much of the work that we have described above leaves 
basic questions unanswered, requiring future observational 
and theoretical work. First, high spatial resolution, from 
microwaves through 7 -rays, is the observational key to un- 
derstanding impulsive-phase transport. What is the spatial 
scale of particle beams? Return currents? To what extent are 
they filamentary? What processes dominate the physics of 
return currents? What physical processes could disrupt 
beams? What role do they play in heating the ambient 
plasma? Does the relative size of hard and soft X-ray sources 
suggest an important energy transport role for microwave 
maser radiation? Second, what is the mass and charge com- 
position of flare-accelerated particles? Do proton beams exist 
in the deka-keV range? Are energetic protons important to 
impulsive-phase momentum transport? Do an energetically 
significant number go undetected because of their inefficiency 
of bremsstrahlung production? Third, where are nonthermal 
particles accelerated within flare loops? What is their degree 
of anisotropy and their energy spectrum? How are they con- 
fined to the corona in those flares that do not show footpoint 
brightening in hard X-rays? How are they thermalized as they 
travel along flare loops? Is their transport dominated by par- 
ticle collisions or by wave generation or by the emission or 
absorption of radiation? Are regions of particle precipita- 
tion identical to regions of explosive coronal and chro- 
mospheric hydrodynamic motion? Will the commonly 
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accepted model of coronal loops as the basic building block 
of flare geometry withstand close scrutiny? Not until the 
answers to these questions are better understood can we claim 
to understand the basic physics of impulsive-phase transport. 
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